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SIGNIFICANT  STRESSES  OF  ARRESTED  SALINE  WEDGES 


PART  I:  INTRODUCTION 

1.  In  1946,  an  experimental  program  on  the  model  laws  of  density 
currents  was  initiated  at  the  National  Bureau  of  Standards  for  the  U.  S. 
Army  Engineer  Waterways  Experiment  Station.  The  experimental  work  con¬ 
tinued  until  1961.  Results  of  various  studies  have  been  reported  with 
the  exception  of  two  studies--one  relating  to  the  interfacial  stresses  of 
the  arrested  saline  wedges  and  the  other,  the  effect  of  the  Richardson's 
number  on  mixing.  The  present  paper  takes  up  the  question  of  stresses. 

2.  Laboratory  studies  are  needed  on  stresses  of  arrested  saline 
wedges  to  provide  guidance  in  the  interpretation  of  physical  models  for 
the  motion  of  saline  waters  into  water  courses  free  of  large  tidal  in¬ 
fluences.  Constructing  the  model  on  the  Freudian  scale,  it  is  necessary 
that  the  relative  densities  be  retained.  In  small  models  the  associated 
interfacial  stresses  are  laminar  in  origin  and  are  dependent  on  Reynolds 
number  of  flow  (Keulegan  1955b).  Since  in  all  probability  the  corre¬ 
sponding  stresses  in  nature  would  be  turbulent  in  character,  transference 
from  the  model  to  the  prototype  must  be  made  with  caution. 

3.  The  present  study  reexamines  further  the  results  discussed  in 
the  report  cited.  Additional  information  is  provided  regarding  the  dis¬ 
tribution  of  velocity  in  a  cross  section  of  the  two  layers,  fresh  and 
saline  waters.  Here,  as  in  the  previous  report,  attention  is  focused  on 
the  mean  values  of  the  stresses  at  the  interface  and  at  the  bottom  of  the 
wedge.  Flow  equations  are  depth-integrated  in  both  layers,  the  fresh 
water  and  the  saline  wedge.  Next,  resulting  equations  are  integrated 
lengthwise,  thereby  yielding  expressions  for  the  mean  stresses.  Here, 
the  derivations  are  repeated  with  more  detail  for  clarity.  Further,  to 
examine  the  validity  of  the  results  involving  the  average  stresses,  the 
principle  of  momentum  is  applied  separately  once  to  the  entire  volume  of 
fresh  water  lying  over  the  wedge  and  next  to  the  volume  of  saline  water 
in  the  wedge. 
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PART  II:  APPARATUS  AND  PROCEDURES 


4.  The  apparatus  consists  of  a  river  forebay,  a  sea,  and  a  channel 
connecting  the  two.  The  horizontal  channel,  400  ft*  long,  2  ft  high, 
and  9  in.  wide,  is  rectangular  in  cross  section.  The  bottom  consists  of 
planed  and  painted  wood.  The  sides  are  of  glass  or  lucite  plates, 
inserted  into  rectangular  steel  frames  which  in  turn  with  the  wooden 
bottom  between  them  are  attached  to  a  steel  base.  The  frames  are  4  ft 
long  and  2  ft  high  and  the  entire  channel  assembly  rests  on  4-ft-high 
concrete  columns. 

5.  Over  the  forebay  a  metering  tank  is  mounted  to  control  the 
discharge  into  the  channel  representing  the  river.  The  bottom  of  the 
tank  consists  of  a  1/8-in. -thick  brass  plate  provided  with  numerous 
circular  orifices  of  different  sizes.  Discharge  into  the  channel  is 
controlled  by  a  chosen  combination  of  orifices,  their  size  and  number, 
and  the  water-surface  elevation  in  the  tank.  The  orifices  are  calibrated 
in  place  by  measuring  the  volume  of  water  collected  in  the  channel 
during  a  known  time  period  with  the  sea  end  of  the  channel  closed.  The 
forebay  is  square  in  shape,  each  side  is  4  ft  long,  and  the  bottom  is 
aligned  with  the  river  channel  bottom.  The  entrance  into  channel  is 
rounded . 

6.  The  sea  is  the  most  critical  part  of  the  entire  apparatus. 

The  10-ft-long,  8-ft-wide,  and  2-ft-deep  rectangular  tank  representing 
the  sea  is  elevated  with  its  horizontal  bottom  aligned  with  the  channel 
bottom.  Despite  the  limited  dimensions,  the  sea  is  required  to  behave 
as  an  unbounded  sea.  It  is  necessary  to  resort  to  a  separate  saline 
water  reservoir  placed  directly  below  to  establish  a  steady  exchange  of 
saline  water.  Return  of  saline  water  from  the  sea  is  governed  by  ele¬ 
vated  crests  of  considerable  length  which  also  maintain  a  selected  sea 
water-surface  elevation.  During  a  test  run  there  is  considerable  mixing 
between  fresh  and  saline  waters  at  the  sea  surface  and  to  maintain  a 


Multiply  feet  by  0.3048  to  convert  to  metres;  multiply  inches  by 
25.4  to  convert  to  millimetres. 
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constant  salinity  of  the  sea,  saline  water  from  the  reservoir  below  is 
constantly  introduced  into  the  sea.  Previous  observations  have  indicated 
that  if  p  +  Ap  be  the  density  of  seawater  and  p  +  Ap'  the  density  of 
fresh  water  at  the  locale  of  efflux  from  the  sea  area,  as  the  result  of 
mixing,  the  ratio  Ap'/Ap  equals  0.34  (Keulegan  1957b).  This  surpris¬ 
ingly  high  value  of  mixing  requires  that  the  density  of  the  saline  water 
in  the  reservoir  be  maintained  at  a  constant  value  through  the  continued 
addition  of  salt  (Figure  1). 


Figure  1.  Plane  view  of  experimental  sea 


7.  Saline  water  densities  are  determined  electrically  by  a 
probe,  the  exploring  end  consisting  of  two  parallel  copper  or  platinum 
wires  2  cm  in  length.  The  salinity  of  an  electrolyte  is  related 
indirectly  to  the  resistance  between  the  two  parallel  wires  immersed 
in  a  solution.  If  the  ratio  of  the  diameter  of  the  wires  to  the  spacing 
between  them  is  small,  the  resistance  is  given  by  (Keulegan  1949): 
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R  =  effective  resistance,  ohins 

T  =  specific  resistance  of  the  electrolyte,  ohms-cm 
L  =  length  of  the  individual  wires,  cm 
s  =  spacing  between  wires,  cm 
d  =  diameter  of  the  wires,  cm 

Since  the  relation  applies  accurately  only  for  wires  of  great  length,  it 

is  used  for  the  purpose  of  designing  the  probe  and  selecting  the  ammeter. 

The  probe  is  calibrated  against  a  standard  solution  of  relative  density 

Ap/p  =0.10  .  The  standard  solution  is  prepared  on  the  assumption  that 

if  V  is  the  volume  of  salt  of  density  p  to  be  dissolved  in  water 
s  s 

of  volume  ,  the  volume  of  the  solution,  all  at  room  temperature,  is 

V  +  V  .  On  this  basis,  the  density  of  the  solution  would  be 
w  s  1 

p  +  Ap  =  (psVs  +  pVw)/(Vg  +  Vw)  (2) 

Lower  density  solutions  are  obtained  by  adding  water  to  standard  solution. 
The  calibration  consists  of  noting  the  ampere  versus  density  when  the 
a-c  voltage  applied  across  the  wires  in  6  volts. 

8.  The  experimental  procedure  for  the  eventual  formation  of  an 
arrested  saline  wedge  is  rather  simple.  In  the  beginning  of  a  test  the 
sea  is  isolated  from  the  river  by  means  of  a  gate.  The  side. gate  at 
one  of  the  channel  walls  about  1  ft  away  from  the  sea  is  left  open 

to  permit  the  lateral  efflux  of  the  river  water.  The  gate  is  manipu¬ 
lated  to  establish  the  desired  depth  of  water  for  the  channel.  With 
this  selection,  the  waters  in  the  channel  and  the  saline  water  of  the 
sea  are  at  the  same  level.  With  this  condition  realized,  the  side  gate 
is  pushed  down  and  the  sea  gate  is  pulled  up  simultaneously.  This  al¬ 
lows  the  entering  of  the  river  current  into  the  sea  area. 

9.  In  the  initial  stages  of  the  movement  the  velocity  of  the  ad¬ 
vancing  saline  front  is  large.  With  distance,  the  velocity  decreases 
gradually.  In  the  initial  stages  the  front  is  rounded  and  also  elevated. 
At  some  critical  distance,  the  roundness  of  the  front  disappears  and  the 
front  assumes  the  appearance  of  a  slowly  moving  wedge.  This  form  of  the 
wedge  is  maintained  even  after  the  saline  wave  is  arrested. 


PART  III:  BASIC  DATA  OF  THE  EXPERIMENTS 


10.  Experiments  were  restricted  to  two  river  depths, 

=  45.5  cm  and  =  23.0  cm  (Figure  2);  these  are  the  depths 


x 


Figure  2.  Notation  diagram,  arrested  saline  wedge 

of  water  immediately  upstream  of  the  tip  of  the  arrested  saline  wedge. 

Since  B,  the  channel  width,  is  22.9  cm,  the  corresponding  depth-width 

waters  are  2  and  1,  respectively.  In  these  tests  the  excess  density  of 

the  sea  saline  waters  over  the  density  of  fresh  water,  Ap  ,  falls  in 

the  range  0.01  to  0.08.  The  river  velocity  just  upstream  of  the 

wedge  is  selected  such  that  n  ,  the  ratio  of  the  depth  of  saline  waters 

of  wedge  at  the  river  mouth  hgl  to  water  depth  ,  is  close  to  1/2. 

The  data  of  observations  regarding  the  river  velocity  ,  the  density 

difference  Ap  ,  the  length  of  arrested  saline  wedge  x^  ,  the  depth  of 

saline  waters  at  the  river  mouth  h  ,  ,  the  fall  of  surface  waters  at 

s  1 

the  river  mouth  AH  ,  and  water  temperature  0  are  shown  in  Tables  1 

and  2  for  the  two  test  series.  Determination  of  the  average  interfacial 

stress  T^  ,  stress  averaged  along  the  entire  length  of  the  interface, 

and  the  average  bottom  stress  T  will  be  evaluated  on  the  basis  of 

o 

the  entries  shown  in  these  two  tables.  The  values  of  Ap  in  Table  1 

are  10  percent  smaller  than  the  original  measurements  and  the  values  of 

h  shown  are  1-1/2  cm  less  than  the  original  observations.  The  reason 
si 

for  these  adjustments  will  be  given  subsequently. 
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PART  IV:  LENGTH  OF  ARRESTED  SALINE  WEDGES 


11.  Expressions  for  the  length  of  the  arrested  saline  wedges  and 
the  depth  of  saline  waters  at  the  river  mouth  in  laboratory  channels  were 
established  previously  (Keulegan  19571)) .  Thus,  a  comparison  with  these 
previous  results  may  indicate  if  the  value  of  some  particular  observa¬ 
tion  of  this  study  requires  adjustment. 

12.  The  length  of  an  arrested  saline  wedge  on  the  basis  of  dimen¬ 
sional  analysis  would  be 


where  V  is  the  river  velocity  and 
r 


(4) 


The  ratio  V  /V.  is  referred  to  as  the  densimetric  Froude  number  and 
r  A 

V^H^/v  ,  as  the  densimetric  Reynolds  number.  When  B/H^  equals  unity, 
the  wedge  length  is 


in  the  densimetric  Reynolds  number  of  range  1  x  10  to  10.0  x  10  .  In 

the  present  study  the  determination  of  length  for  the  case  Hq  =  23.0  cin 

and  B  =  22.9  cm  on  the  basis  of  values  shown  in  Table  2  shows  agreement 

with  the  above  as  the  mean  value  of  the  coefficient  is  also  0.19.  When 

U 

B/H  equals  2  in  the  densimetric  Reynolds  number  of  range  4.0  x  10 
°  4 

to  40.0  x  10  the  length  is  given 
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When  the  determinations  of  length  for  the  case  H  =  45 . 5  cm  and 

o 

8  =  22.9  cm  were  made  on  the  basis  of  the  recorded  values  of  Ap/p  , 

which  were  10  percent  higher  than  those  shown  in  Table  1,  the  results 

were  not  in  agreement  with  Equation  6;  the  mean  value  of  the  coefficient 

was  0.10.  Repeating  the  determination  with  Ap/p  values  as  shown 

in  Table  1,  the  mean  value  of  the  coefficient  is  0.12  in  agreement  with 

Equation  6.  Due  to  error  of  probe  calibration,  sea  densities  are  not 

correctly  determined.  The  error  does  not  affect  the  evaluation  of  the 

interfacial  mean  stress  T  ,  but  it  would  affect  the  bottom  mean  stress 

s 

T^  ,  using  the  formulae  subsequently  established. 

13.  Previous  studies  (Keulegan  1957b)  have  indicated  that  the 

depth  of  seawater  at  the  river  mouth  expressed  as  ratio  h^/H^ 

practically  independent  of  the  density  of  seawater,  the  width-depth 

ratio  of  the  channel,  and  the  densimetric  Reynolds  number;  it  is  merely  a 

function  of  the  densimetric  Froude  number  V  /V,  as  shown  in  Table  3. 

r  A 

The  entries  represent  averages  from  tests  with  different  channels  and 
different  width-depth  ratios  (Keulegan  1957b). 

14.  As  first  shown  by  Shijf  and  Schoenfeld  (1953),  the  flow  of 
fresh  water  at  the  river  mouth  is  critical, 


U 


2 

1 


Ap 

P 


ghwl 


1 


(7) 


where  h  ,  is  the  depth  of  fresh  water  at  the  river  mouth  and  U,  is 
wl  r  1 

the  velocity  of  freshwater  current  (Keulegan  1955a).  Assuming  that  the 
fall  of  surface  water,  AH  ,  is  negligible  one  has  h  ,  =  H  -  h  ,  and 
thus 


(8) 


The  entries  in  Table  3  are  in  good  agreement  with  this  expression  for 
V  /V  close  to  0.5  only.  For  smaller  values  of  V^/V^  >  t-heory  yields 
depth  values  greater  than  those  given  in  this  table. 
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15.  In  the  derivation  of  the  theoretical  expression,  the  theory 
assumes  tacitly  that  the  distribution  of  pressure  is  hydrostatic.  On 
the  other  hand,  when  the  densimetric  Froude  number  is  small  the  inter¬ 
face  in  the  area  of  river  mouth  would  be  curved,  a  condition  wherein  the 
pressure  would  be  no  longer  hydrostatic.  Thus,  it  is  desirable  to  obtain 
the  depth  hg^  from  a  relation  based  on  observation. 

16.  In  this  study  the  position  of  the  interface  was  obtained  visu¬ 
ally.  To  enhance  a  stronger  contrast  between  the  fresh  and  saline 
waters,  the  water  of  the  sea  was  colored  by  introducing  a  chrome  com¬ 
pound.  This  is  ideal  if  there  is  no  mixing  at  the  interface.  Un¬ 
fortunately,  if  agitation  is  present  small  traces  of  this  material  may 
affect  a  deep  coloring  and  the  true  position  of  the  interface  may  not 

be  inferred  correctly.  This  was  the  case  with  the  tests  carried  with 

water  depth  H  =  45.5  cm  and  thus  the  inferred  values  of  h  ,  are  not 
r  o  si 

reliable.  Indeed,  the  ratio  Hg^/H  when  computed  from  the  inferred 
values  are  at  variance  with  the  entries  of  Table  3.  However,  agreement 
would  be  obtained  only  after  subtracting  1.5  cm  from  the  inferred  h  .  . 

S  X 

These  adjusted  values  of  h  .  are  shown  in  Table  1.  A  similar  adjust- 

S  X 

ment  is  not  necessary  for  the  tests  carried  out  with  water  depth 

H  =  23.0  cm. 
o 
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PART  V:  SHAPE  OF  WEDGES  AND  FAIX  OF  SURFACE  WATERS 


17.  Information  on  the  shape  of  the  arrested  saline  wedge  and  the 
fall  of  surface  waters  over  the  wedge  plays  an  important  part  in  de¬ 
termining  the  average  boundary  stress  at  the  interface  and  at  the  bottom 
from  formulae. 

18.  Previous  studies  have  shown  that  within  a  definite  range  of 
densimetric  Froude  number  the  shapes  are  affine  to  each  other  irrespec¬ 
tive  of  river  velocities,  the  density  of  seawater,  and  channel  width 
(Keulegan  1957b).  Denoting  by  hg  the  depth  of  saline  waters  in  the 
wedge  at  a  point  of  distance  x  from  the  tip  of  the  wedge 


where  denotes  the  length  of  the  wedge. 

Introducing  the  relative  distance 

t  =  —  (101 

x 

o 

the  affine  shape  is 

~  =  f(0  Ol) 

si 

This  functional  relation  is  shown  in  Table  4. 

19.  For  the  present  study  the  observation  of  shape  of  arrested 

saline  wedge,  selecting  the  river  water  depth  of  ratio  of  Hq  =  23  cm, 

is  restricted  to  tests  with  a  common  densimetric  Froude  number 

V  /V.  =  0.40  and  five  different  seawater  densities.  Tests  with  a 
r  A 

given  seawater  density  are  repeated  and  the  averages  of  hg/hsj  are 
sought;  these  are  shown  in  Table  5.  Comparison  between  the  entries  of 
different  columns  indicates  that  for  a  constant  densimetric  Froude 
number  the  affine  shape  is  independent  of  the  river  velocity  or  the 
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20.  To  measure  the  fall  of  surface  waters,  along  the  60-m  length 
of  the  experimental  channel,  12  glass  manometer  tubes  were  mounted 
vertically  against  the  glass  sidewalls.  Next  to  the  tubes,  paper  scales, 
(reading  in  centimetres)  were  glued  to  the  glass.  The  meniscus  positions 
were  read  by  a  telescope.  Displacement  of  the  menisci  was  readily 
estimated  within  a  tenth  of  a  millimetre.  To  establish  a  common  zero 
for  the  scales,  the  exit  gates  of  the  channel  were  closed  and  the  channel 
was  filled  with  fresh  water  to  a  designed  depth  Hq  ,  23.0  cm  or  45.5  cm. 
The  readings  of  the  menisci  in  all  the  manometer  tubes  were  taken  after 
all  the  surface  oscillations  had  ceased.  An  example  of  the  measurement 
of  the  fall  of  the  freshwater  free  surface  is  given  in  Figure  4.  The 


Figure  4.  An  example  of  observed  values  of 
fall  of  water  surface 


part  of  the  curve  corresponding  to  x  negative  represents  the  fall 
upstream  of  the  tip  of  the  arrested  wedge.  The  slope  here  is  in  accord 
with  the  relation 
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(13) 


Ah  _  .  _~o_ 

Ax  "A  2gR, 


where  is  the  hydraulic  radius  and  \  the  coefficient  of  resistance 


0.0559 


The  part  of  the  curve  corresponding  to  x  positive  represents  the 
surface  fall  of  fresh  water  lying  over  the  arrested  saline  wedge. 

Denoting  the  fall  of  this  surface  at  the  river  mouth  by  AH  and  that  at 
a  point  £  =  x/x^  by  Ah  ,  the  ratio  Ah/AH  will  be  referred  to  as  the 
relative  fall  of  surface  waters.  This  quantity  as  a  function  of  £  is% 
determined  from  curves  similar  to  the  one  shown  in  Figure  4.  Determina¬ 
tions  made  with  varying  seawater  densities  and  for  a  depth  of  water  = 
23.0  cm  are  shown  in  Table  6  and  those  corresponding  to  Hq  =  45.5  cm 
in  Table  7.  It  appears  that  within  the  range  of  the  densimetric  Froude 
number  used  the  ratio  is  independent  of  this  number.  The  average  values 
of  relative  fall  from  Tables  6  and  7  are  transferred  to  Figure  5  and  the 
curve  through  the  points  yields  the  relation 


—  =  £-bsin7t£  +  csin  2 


where  b  =  0.17  and  c  =  0.03  . 

21.  In  the  case  of  an  arrested  saline  wedge  with  the  fresh  water 
above  it  in  motion,  it  would  be  easier  to  evaluate  the  average  inter¬ 
facial  stress  by  directly  considering  the  hydrodynamics  of  the  freshwater 
layer.  For  this  determination  it  is  important  that  both  AH  and  Ah/AH 
are  correctly  observed  and  evaluated. 
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PART  VI: 


VELOCITY  PATTERNS 


22.  The  determination  of  the  mean  stress  T  of  the  interface 

s 

and  the  mean  bottom  stress  T^  on  the  basis  of  the  dynamics  of  the 
arrested  saline  wedge  requires  that  the  pattern  of  the  velocities  in 
this  area  is  accurately  known.  The  problem  becomes  more  complicated  if 
there  is  a  flow  due  to  mixing  across  the  interface. 

23.  Examination  of  velocities  in  the  arrested  wedge  and  in  the 
part  of  fresh  water  lying  over  the  wedge  was  restricted  to  tests  with  a 
river  depth  of  23.0  cm.  Interfacial  velocities,  surface  water  veloc¬ 
ities,  and  velocities  of  fresh  water  and  saline  water  in  a  cross  section 
were  measured  separately. 

24.  For  the  observation  of  the  interface  velocities  a  small 
globule  of  butyl  phthalate  dissolved  in  xylene,  with  a  density  slightly 
less  than  the  density  of  the  seawater,  is  injected  by  means  of  an  eye¬ 
dropper  into  the  fresh  water  a  few  centimetres  below  the  surface.  After 
the  descent  of  the  drop  is  completed  and  it  reaches  the  interface, 
subsequent  motion  of  the  drop  from  one  locality  to  another  is  timpd. 

Let  IF  be  this  interface  velocity  observed  at  the  point  x  and  U  be 
the  mean  velocity  of  fresh  water  in  the  cross  section  through  x  . 
Forming  the  ratio  IF/U  ,  observed  values  from  the  tests  corresponding 
to  various  seawater  densities  are  shown  in  Table  8  as  function  of  £  . 
For  a  given  x/x^  the  ratios  IF/U  reveal  small  differences  from  one 
column  to  the  other.  These  differences  may  be  ignored.  An  interface 
with  a  sharp  discontinuity  in  the  densities  does  not  exist.  Actually, 
there  is  a  gradation  of  densities  between  the  two  liquids  and  one  takes 
the  surface  where  the  density  equals  1  +  1/2  Ap  as  the  effective  inter¬ 
face.  The  density  of  the  butyl  phthalate  globules  cannot  be  accurately 
controlled  and  in  different  tests  the  globules  are  differently  placed 
with  respect  to  the  interface.  It  is  thus  sufficient  to  consider  the 
average  values  entered  in  the  last  column  of  Table  S  The  plotting  in 
Figure  6  is  on  this  basis.  Starting  from  zero,  the  value  of  LF/U 
increases  very  rapidly  and  then  slowly  reaches  asymptotically  the  limit 
U./U  =  0.53  .  The  surface  velocity  is  determined  by  timing  the 
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Figure  6.  Interface  and  surface  velocities,  Fq  =  0.40 

motion  of  small  paraffin  particles  riding  on  the  surface  waters.  Also, 

in  this  case  the  ratio  U  /U  is  formed  and  the  values  as  a  function  of 

s 

£  from  the  tests  using  different  seawater  densities  are  entered  in 

Table  9.  Here  also  U  /U  values  reveal  small  differences  from  one 

s 

column  to  the  other.  The  average  values  are  entered  in  the  last  column 
and  the  plotting  in  Figure  6  is  on  this  basis. 


25.  The  device  to  determine  the  distribution  of  the  velocities  in 
a  flow  cross  section  consisted  of  a  thin  phosphor  bronze  strip  about 
1/4  in.  wide  and  8-1/2  in.  long  with  one  of  the  ends  of  the  strip 
soldered  to  a  rectangular  rod.  Touching  one  of  the  sidewalls  of  the 


experimental  channel,  the  rod  was  held  vertical  with  the  strip  facing 
the  current  normally  in  a  horizontal  plane.  The  deflection  6  of  the 
free  end  of  the  strip  is  an  indication  of  the  velocities  in  the  plane  of 
the  strip.  An  argument  of  dimensional  analysis  yields 


(16) 


where 


6  =  deflection  of  the  free  end  of  the  strip 
SL  and  t  =  length  and  thickness,  respectively,  of  the  strip 
p  =  density  of  liquid 

u  =  root  mean  square  of  the  velocities  in  the  strip  plane 
E  =  Young's  modulus  of  the  strip  natural 
As  the  deformation  of  the  strip  is  flexural  in  the  main 


=  f 


[*  (Ol 


(17) 


For  a  given  strip  E  ,  t  and  iL  are  constants.  Then, 


p1/2u  =  f (6)  (18) 

where  the  quantities  are  measured  in  standard  units,  cgs .  Although 
the  relation  may  be  established  using  Kirchhoff's  theory  of  slender 
wires,  it  is  more  practical  to  resort  to  calibration.  For  this  purpose 
the  strip  is  held  at  middepth  of  a  current  in  the  experimental  channel 
of  discharge  Q  and  the  end  deflection  6  is  noted.  The  average 
velocity  u  is  deduced  from  the  discharge  and  the  strip  elevation  from 
the  channel  bottom  using  Blasius'  velocity  relation.  It  is  assumed 
that  the  velocities  of  all  the  points  in  a  cross  section  having  the 
same  common  elevation  are  equal.  This  introduces  a  small  calibration 
error  since  the  velocities  at  the  points  of  strip  level  are  not  constant 
This,  however,  will  not  severely  affect  the  value  of  u/U  ,  where  u  is 
the  root-mean-square  value  of  the  velocities  at  points  of  a  common  eleva 
tion  from  the  bottom  and  U  is  the  average  velocity  of  the  freshwater 


current  in  the  same  cross  section,  all  the  velocities  being  baser!  on  a 
chosen  calibration. 

26.  Using  five  different  seawater  densities  the  examination  of 

velocities  was  confined  to  a  central  area  of  x/x  =  0.75  where  the 

o 

depth  of  saline  water  was  7.0  cm  and  that  of  the  fresh  water  above, 

16.0  cm.  All  the  arrested  wedges  realized,  corresponding  to  a  densi- 

metric  Froude  number  F  =  0.40  ,  were  nearly  of  a  common  length  x  = 

o  J  o 

2000  cm. 

27.  Relative  velocities  u/U  of  the  saline  wedge  area  are  given 
in  Table  10  and  those  of  fresh  water  over  the  wedge  in  Table  11.  The 
entered  values  were  read  from  individual  curves  drawn  smoothly  through 
points  of  observation.  Actually,  there  was  some  scatter  in  the  observed 
values.  The  small  variation  in  the  entry  values  from  column  to  column 
suggests  that  the  pattern  of  velocities  in  a  cross  section  is  hardly 
affected  by  the  seawater  densities,  all  the  arrested  wedges  corresponding 
to  a  common  densimetric  Froude  number.  Figure  7  was  prepared  using  the 
mean  values  shown  in  the  last  column  of  Tables  10  and  11. 

28.  Flow  in  the  tip  area  of  the  arrested  saline  wedge  with  £ 
less  than  0.1  may  be  supposed  to  be  irrotational ,  the  same  as  in  the 
case  where  the  boundary  is  composed  of  two  rigid  walls  inclined  at  an 
angle  less  than  90  deg.  In  the  remaining  part  of  the  wedge  if  there  is 
mixing  at  the  interface,  this  will  influence  the  velocities  in  the  lower 
areas  of  the  arrested  wedge  where  the  motion  is  away  from  the  sea. 

Using  the  methods  of  an  earlier  investigation  on  mixing  in  arrested 
saline  wedges  (Beta  1957),  it  may  be  estimated  that  the  total  quantity 
of  saline  water  traversing  the  interface  into  fresh  water  due  to  mixing 
is  about  8  percent  of  the  saline  water  inflow  into  the  wedge  from  the 
sea.  The  effect  of  this  may  be  ignored.  One  may  then  assume  that  the 
velocity  patterns  from  one  cross  section  to  another  are  similar  to  each 
other.  It  was  already  seen  that  the  interfacial  velocity  U^/U  has 
practically  the  same  value  along  almost  the  entire  length  of  the  inter¬ 
face  except  in  the  area  near  the  tip. 
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FART  VII:  KYDRiWNAMI  CS  OF  FRESHWATER  AND  INTF.RFAC7AL  STRESS 


29.  Let  the  origin  of  the  rectangular  axis  system  be  placed  at 
the  tip  of  the  arrested  saline  wedge  and  at  the  midpoint  of  the  channel 
bottom  width.  Let  the  axis  of  x  be  drawn  horizontally  in  the  di¬ 
rection  of  motion  of  the  fresh  water  and  z  vertically  upward.  Com¬ 
ponents  of  the  velocities  along  the  axes  x  ,  y  ,  and  z  are  u  ,  v  , 
and  w  ,  respectively.  Elevation  of  the  water  surface  measured  from 
the  horizontal  channel  bottom  may  be  denoted  by  h  and  that  of  the 
interface  by  h  .  Let  the  width  of  the  channel  be  B  ,  B  =  2b 
(Figure  8) . 


Figure  8.  Notation  diagram,  analysis  of  wedge  stresses 

30.  As  the  length  of  the  arrested  saline  wedge  is  large  in  com¬ 
parison  with  the  maximum  depth  h  ^  of  the  wedge  at  the  river  mouth, 
the  vertical  component  of  the  fluid  acceleration  will  be  neglected  and 
the  pressure  will  be  evaluated  hydrostatically.  Under  this  assumption, 
using  the  notation  of  Lamb,  the  equations  of  motion  are 


3u  9u  1  3p  1  /  ^zx  ^yx  \ 

U  3x  3z  ~  '  p  3x  p  y  3z  3y  J 


(19) 
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+b  h 

J  j  [a!  (u2)  +  b  (uw)  dz  <ly] 

-b  hs 
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where  the  subscript  "s"  refers  to  the  surface  and  "i"  to  the  interface. 
As  a  particle  on  the  water  surface  stays  oh  this  surface  and  a  particle  1 
on  the  interface  stays  on  the  interface  in  the  absence  of  mixing 


w 

s 


u 


dh 
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dh 

_ s 
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Then,  the  final  form  of  the  left-hand  member  of  Equation  21  is 


* 


a  a  ^  f  2  dh 
dz  dy  +  u  — 
s  dx 


2 
u  . 
i 


(22  1 


One  now  introduces  a  ,  the  Boussinesq  coefficient  of  velocity  distribu¬ 
tion,  defined  by 


+b 


where  U  is  the  average  value  of  u 
Differentiating  the  two  sides  of  Llii 


h 

u2  dx  dz  (23  ) 

h 

s 

in  a  cross  section  A  =  B(h  -  h  ) 

s 

last  equation  with  respect  to  x  , 
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in  accordance  with  thp  l.eihnit/  rule, 


+b  h 

d  ,  2  .  ,  /  /  a  2.  /  2  dh  2  dh  \  D 

-  (aU  A)  =  J  J  8x  u  dz  dy  +  -  -  u.  ^  j  B 

-h  h 

s 

Hence  after  comparing  with  Equation  22,  the  left-hand  member  of  Equa¬ 
tion  21  is 

h  +b 

J  f  (uk  +  uiz)dzdy  =  k  (m2a)  (24 ’ 

h  -b 
s 

A  slight  variation  of  a  with  respect  to  x  may  be  ignored  and  as  UA 
is  constant 

55  Ua)  =  UA  55  (au)  -a2At2  =  l  €  ih-  hs)B  (25) 

Thus  finally 

+b  h 

/[  (  9u  A  3u  \  _  a  dU2  ,  . 

J  [ud^  +  w  Vz)dz  dy  -  2  d5  (h  '  hs)B 

-b  h 

s 

This  is  an  evaluation  of  the  integral  on  the  left-hand  side  of  Equa¬ 
tion  21.  From  Equation  20 

P  =  Pa  +  Pg(h  -  z)  h  -  hs 
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where  p  is  the  atmospheric  pressure.  Accordingly,  the  first  integral 

3  ” 

on  the  right-hand  side  of  Equation  21  reduces  to 
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(27) 


In  the  absence  of  air  currents,  the  stress  at  the  water  surface  vanishes. 

Let  the  traction  of  the  fresh  water  on  the  saline  water  be  denoted  by 

t  .  Thus 
s 
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zx 
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z  =  h 
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and  the  second  integral  on  the  right-hand  side  of  Equation  21  reduces  to 


I 

P 


(28) 


where  Tg  is  the  average  value  of  the  interfacial  stress  along  the  chan¬ 
nel  width.  It  is  a  positive  quantity.  Denoting  the  resistive  force  of 
the  vertical  wall  by  ,  a  positive  quantity, 
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the  third  integral  on  the  right-hand  side  of  Equation  21  reduces  to 

b  h 

ill 

-b  h 


where  t  is  the  average  value  of  the  wall  shear  across  the  span 
w 

h  -  hs  .  In  view  of  Equations  26,  27,  28,  and  29,  the  original  Equa¬ 
tion  21  reduces  to 
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(30) 


This  is  the  expression  for  the  local  stress  in  terms  of  the  surface 

slope,  the  surface  elevation,  the  saline  wedge  height,  and  the  current 

mean  velocity.  All  these  refer  to  the  situation  at  the  point  x  . 

However,  owing  to  the  difficulties  in  measuring  the  surface  fall  rate 

dh/dx  accurately,  it  will  prove  to  be  more  serviceable  to  consider 

the  average  T  of  the  interface  stress,  averaged  along  the  entire 
s 

length  of  the  saline  wedge. 
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(31) 


and  in  view  of  Equation  30 
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Since 

U  H  =  U, (H  -  AH  -  h  ,)  =  U(h  -  h  )  (33) 

o  1  o  si  s 

where  is  the  velocity  of  fresh  water  at  the  river  mouth  and  AH 

toe  fall  of  surface,  the  last  term  in  Equation  32  changes  to 
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(H  -  AH  -  h  .)  -  U2H 
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Also,  neglecting 
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(34) 
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For  brevity  introduce  T  , 

wl 

the  fresh  water 


the  total  frictional  force  of  the  walls  on 


T 


wl 


(37 ) 


In  view  of  the  last  four  relations,  Equation  32  changes  to 


32.  One  may  also  obtain  this  relation  by  applying  the  principle 
of  momentum  to  the  body  of  water  found  at  the  instant  t  over  the  wedge 
with  its  boundary  ABDE  as  shown  in  Figure  9.  Let  dM/dt  denote  the 
rate  of  change  of  momentum  of  this  water  in  the  direction  x  positive. 


Figure  9.  Notation  diagram  relating  to  momentum  of  fresh  water 
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The  forces  producing  this  change  are  the  pressure  forces  and  ?  ^ 

acting  on  the  upstream  and  the  downstream  forces  AB  and  ED,  and  the  pres¬ 
sure  force  P.  on  the  interface.  The  tangential  traction  forces  taking 
part  in  this  change  of  momentum  are  T^^  ,  the  retracting  force  from  the 
two  sidewalls,  and  TgBxo  ,  the  retracting  force  due  to  the  saline  wedge. 


tt  =  P,  -  P.  -  P.  -  T  ,  -  T  Bx 
fit  1  2  1  wl  so 


The  particles  of  water  which  at  time  t  are  in  the  planes  AB  and  DE  move 
into  the  surfaces  A'B  and  D'E'  at  time  t  +  fit  ,  respectively.  Let  Mq 
be  the  momentum  of  the  liquid  contained  in  A'BDE;  Mj  of  that  contained 
in  ABA';  and  of  that  contained  in  EDD'E'.  On  this  basis  the  momen¬ 

tum  of  liquid  under  consideration  is  M  +  M,  at  the  instant  t  and 
Mo  +  at  the  instant  t  +  fit  .  Thus  the  rate  of  change  of  momentum  is 


^  =  (M2  -  Mj)fit 
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where  Uq  is  the  river  current  average  velocity  at  the  tip  of  the  saline 
wedge  and  the  average  velocity  of  fresh  water  at  the  river  mouth. 
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Introducing  these  in  Equation  40  and  assuming  =  a 


u  =  p“  [ui2<"„  -  M  -  »5l>  -  "X  B] 


(41) 


Ignoring  the  pressure  of  the  atmosphere,  the  end  pressure  forces  are 


pi  =  p*  r  8 


and 


Thus,  ignoring  AH 
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P2  =  Pg  -2 - y— *3-  B 
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The  total  pressure  force  from  the  interface  is 

h  i 

sl 

P.  =  pgB  /  (H  -  Ah  -  h  )  dh 


/ 


O  OSS 


or 


Now 


_d 

dx 


dh 

(Ahhs)  =  ih  jji  t  hs 


dPh 

dx 


30 


and  as 


finally 


From  Equations  42  and  43 


(43) 


Substituting  from  Equations  41  and  44  in  Equation  39  and  dividing  the  re¬ 
sulting  equation  by  pBxQ 


This  last  equation  is  identical  with  Equation  $8,  the  alternate  form  of 
Equation  32  which  was  derived  by  integrating  vertically  the  Eulerian 
flow,  Equation  19. 
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PART  VIII:  HYDRODYNAMICS  OF  ARRESTED  SALINE  WEDGE 
AND  AVERAGE  BOTTOM  STRESS 


33.  Multiplying  the  terms  in  Equation  19  by  dz  dy  and  inte¬ 
grating  over  the  cross  section  of  the  wedge  at  x  , 


1_ 

P’ 


+b  h 


0£ 

3x 


dz  dy 


+b  h 


b  o 


dz  dy  +  ^r 
P 


+b  h 


-b  o 


dz  dy 


(46) 


where  p'  is  the  density  of  saline  waters  of  the  edge.  Since  in  the 
normal  cross  section  of  the  arrested  saline  wedge  the  mean  velocity  U 
practically  vanishes,  one  now  introduces,  in  establishing  an  expression 
similar  to  Equation  23,  the  interfacial  velocity  IL  and  defines  a  new 
coefficient  of  velocity  p  as 


PU*A 


+b  h 


b  o 


dz  dy 


(47) 


Using  the  same  argument  that  was  applied  to  the  freshwater  part  it  may 
be  shown  that 


(48) 
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Assuming  that  the  pressure  is  hydrostatic  also  in  the  wedge 


p  =  p  +  Pg(h  -  h  )  +  g(p  +  Ap)(h  -  z)  ,  z  <  h 

3  S  S  S 

where  Ap  is  the  excess  of  density  of  the  saline  water  of  the  wedge 

over  that  of  fresh  water;  p'  =  p  +  Ap  .  Differentiating  with  respect 
to  x 


= 

3x 


3h 
,P  3x 


+  Ap 


3h 
_ s 

3x  , 


and  setting  in  the  first  integral  of  the  right-hand  side  of  Equation  46 


h  +b 


1  3p 
p’  3x 


a  a  (  ^ 

dz  dy = 


a  dh  \ 

__s ) 

P  3x  / 


h  B 
s 
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since  p'/P  is  practically  unity.  Denoting  the  frictional  stress  of  the 
bottom  by  r 

o 


P 


zx 


t 


s 


z  =  h 

s 


P 


zx 


-T 

o 


z  =  0 


the  quantities  x  and  Tq  being  regarded  as  positive.  The  second 
integral  on  the  right-hand  side  of  Equation  42  now  becomes,  neglecting 
again  the  difference  between  p  and  p'  , 

+b  hg  +b 

K  I  I  -Ir*  dz  dy  =  -  /  (T  +  T)dy  =  +  -(i  +  l)B  (SO) 

p  II  dz  PI  s  o  3  p  \  s  o } 

-b  o  -b 

In  regard  to  the  sidewall  effect 
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P  =  "T 
yx  w 


y  =  h 


P 


yx 


+  T 

w 


y  =  -b 


Here  is  positive  when  the  flow  in  the  neighborhood  of  the  wall  is  in 

the  direction  of  the  main  flow  of  fresh  water  above  and  negative  when  in 
the  opposite  direction.  The  sign  of  I  changes  from  the  positive  to  the 
negative  in  moving  from  the  interface  downward.  Now, 


+b  h 


b  o 


dz  dy 


T  dz  =  -  -  iwh 
w  p  s 


(51) 


Substituting  the  expression  from  Equations  48,  49,  50,  and  51  in  Equa¬ 
tion  46,  dividing  by  B  ,  and  transferring  terms 


+  !o  =  /|h  +  Ag  ^s\  h  +  2  -  j>  +  d 

p  p  \8x  p  dx  /  s  p  w  B 


k  (Khs) 


(52) 


Multiplying  the  terms  in  this  equation  by  dx  ,  integrating  between  the 
limits  x  =  0  and  x  =  xq  ,  and  writing 


x 

r  - 

r  x  =  /  t 

o  o  /  0 


dx 


(53) 


T  T 

s  .  o  _ 

—  +  —  =  g 

P  P 


x 

/' 


dh  .  dx  1  Ap 

dx  s  x  2  8  p 

o 


h2l  2f° 

Xo  Pj  w 


s  dx 
B  x 


pu2,  ^ 

1  1  X 

o 


(54) 


This  last  relation  can  be  obtained  also  by  applying  the  momentum  principle 
to  the  entire  body  of  the  salt  wedge  at  the  instant  t  and  contained  in 
the  boundary  ADCB  as  shown  in  Figure  10.  Let  the  rate  of  change  of 
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Figure  10.  Notation  diagram  relating  to  momentum  of  salt  wedge 


momentum  of  this  body  of  liquid  in  the  direction  of  x  positive  be  de¬ 
noted  by 

6M 

6t 


The  forces  producing  this  change  are  the  total  pressure  forces  P. 
acting  on  the  interface  and  acting  on  the  limiting  face  of  the 


saline  wedge  with  the  depth  h 


s  1 


The  tractive  forces  also  taking  part 


in  this  are  the  interfacial  total  stress  T  Bx  and  the  retarding 

s  o 


force  of  the  sidewalls 
Thus 


w2 


~  =  P.  -  P  +  T  Bx  +  T  Bx  -  T 
ot  1  2  so  oo  w2 


(55  ) 


In  the  notation  previously  adopted 


T 


w2 


o 


dx 


(56) 
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The  particles  of  saline  water  which  at  time  t  are  in  the  planes  AO  and 

OC ,  move  into  the  surfaces  ADO  and  OC ' ,  respectively ,  at  time  t  +  6t  . 

Let  be  the  momentum  of  liquid  in  ADOCB,  of  that  contained  in 

ADO,  and  of  that  contained  in  OC ’ C .  On  this  basis  the  momentum  of 

the  saline  bodv  of  water  under  consideration  is  M  +  M,  at  time  t 

o  1 

and  M  +  M„  at  the  later  instant  t  +  6t  .  Thus  the  rate  of  change  of 
o  2 

momentum  is 


Now, 


<5N 

6t 


(M2  -  Mj) 
6t 


(57  ) 


6  +b 


o  -b 


dz  dy  6t 


and 


dz  dy  6t 


as  Uj  is  negative  for  z  smaller  than  6  ,  and  positive  for  z 
greater  than  6  .  Substituting  in  Equation  53 


h  +b 

«  J  f  2  „  „ 

6t  =  BP/  J  "]  dx  d>’ 

o  -b 

and  in  agreement  with  Equation  47,  A  =  h^B 


(58) 


fiM  -  RII2h  R 

6t  ~  PVslB 


(59) 
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From  Equation  39,  rewriting 


where 


P2  =  pg  (Hq  -  AH  -  hgl)  +  g(p  +  Ap)  hgl  -  g (p  +  Ap)  z 


The  integration  yields 


P 


2 


PgB 


-  AHh  , 
si 


'll  ,4e  ii) 

2  P  2  j 


(61.) 


Substituting  in  Equation  51  from  Equations  59,  60,  and  61  and  dividing  by 
pBx^  ,  results  in 


T 

s 

P 


/  h 

dAh  dx  +  g  Ag  s 1 

s  dx  x  2  p  x 
o  o 


(62) 


which  agrees  with  Equation  54,  since  h  =  H 


Ah  . 
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PART  IX:  DIMENSIONLESS  FORM  OF  THE  STRESS  EQUATIONS 


34.  In  evaluating  the  average  interface  and  bottom  stresses  on 
the  basis  of  Equations  32  and  50  using  experimental  data,  it  will  be  use¬ 
ful  first  to  express  these  equations  in  a  dimensionless  form.  It  is 
further  required  that  the  average  resistive  stress  2r  on  the  channel 
walls,  from  fresh  water  and  saline  water,  be  ascertained  analytically 
since  these  stresses  are  not  amenable  to  direct  measurements. 

35.  In  laboratory  experiments  the  Reynolds  number  of  flow  ac¬ 
counted  is  generally  small  and  hence  the  resistance  of  the  sidewalls 
may  be  determined  by  the  Blasius  relation 

1/7 

,  a  =  7.64  (63) 

1/2  - 

where  u.v  =  (t/p)  and  u  is  the  mean  velocity  at  points  having  a 

constant  distance  z'  from  the  interface.  From  this 


and 


0.0284 


where  U  is  the  average  velocity  of  the  fresh  water  in  a  section. 
Then  for  the  average  shear  on  the  wall 


1 


o 
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z'  is  distance  measured  from  the  interface  and  h  is  the  depth  of 
fresh  water.  First  affecting  the  integration  on  the  basis  of  data  shown 
in  Table  11,  one  then  has,  after  introducing  ,  the  river  mean 

velocity , 


2x 


(64) 


where 


A 

o 


0.051 


(65) 


This  estimate  is  only  provisional  as  it  assigns  a  larger  value  to  the 
resistance.  The  derivation  ignores  the  effects  of  turbulence  of  the 
free  surface  and  of  the  interface.  At  such  points  one  expects  lower 
stresses  than  assumed.  In  the  saline  water  area  the  frictional  effect 
of  the  wall  is  not  very  critical.  Since  in  one  part  motion  is  directed 
toward  the  sea  and  in  the  other  part  away  from  the  sea,  resistance 
effects  are  somewhat  neutralized.  As  the  application  of  the  above  method 
to  this  case  is  less  likely  to  be  valid,  it  is  better  to  ignore  it  for 
the  present. 

2 

36.  Dividing  Equation  32  by  ,  the  result  is 


where 


gH 


o  AH 


o 


X 

o 


o  B 


(66  ) 


d(dh/AH) 

dC 


,  n 


H 
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h  +  AH 
s  1 


si 


2  (H  -  AH 
o 


h  .) 
s  1 


H  -  h  , 
o  s  1 
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and 


7K 


-7/4  / h  -  h 


H 


d£  = 


] 

/ 


v-3/4 

'•vl 

s  ] 


(671 


Utilizing  the  experimental  values  of  hg/hsl  and  Ah/A/H  expressed  in 
terms  of  £  ,  see  Equations  12  and  15,  and  affecting  the  integration  and 
remembering  that  n  =  h  /H  >  the  above  given  multipliers  are 


Ij  =  1  -  0 . 59n 


(68) 


I2  =  n/(l  -  n) 


(69 ) 


I3  =  1  +  0 . 37n  +  0 . 20n 


(70J 


37.  Dividing  Equation  54  by  U'  ,  ignoring  the  small  value  term 
from  the  sidewall  fraction,  the  result  is 


T  T  gH  AU  gH  .  h  ,  H 

o  s  o  AH  ,  jn6oAp  sllT  o 

„,,2  -,,2  ..2  x  4  2  ,,2  p  x  5  x 

pU  pU  U  o  U  o  o 

o  o  o  o 


(711 


where 


1 

/ 


T  _  ,  d(Ah/AH)  s  „„ 

1,  -  /  - jp - h~  -  °'59 

si 


(72) 


'5  =  ^V/  (.  -n,2 


The  definition  of  the  Boussinesq  velocity  coefficients  a  is  given  by 
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Equation  23  and  that  of  p  by  Equation  47.  These  coefficients  may  be 
evaluated  through  the  entries  of  the  last  columns  in  Tables  10  and  11, 
respectively.  It  will  be  remembered  that  these  tables  give  the  root- 
mean-square  values  of  velocity  u  of  points  of  as  common  distance  z' 
from  the  interface  in  a  given  cross  section.  This  fact  is  due  to  the 
manner  in  which  the  determinations  are  made.  The  hydrodynamic  forces 
acting  on  the  ribbon  and  causing  it  to  deflect  are  proportional  in  the 
square  of  the  velocity  of  the  particles  striking  the  ribbon.  With  this 
interpretation 

a  =  1.012  (73) 

and 

P  =  0.14  (74) 

since  U.,/U,  =  0.53  and  8  =  0.14  the  last  multiplier  reduces  to 
xl  1  r 

I  =  —  •03-^n  (751 

(1  -  n)2 

Accordingly,  after  neglecting  the  term  involving  \  ,  representing  the 
effect  of  wall  friction  the  final  form  of  the  equation  to  evaluate  the 
average  bottom  shear  is 


T 

o 


+ 


^-0.59 


AH 

X 

o 


+  I  _si)  + 
2  p  x„  / 


0 . 033n 
(1  -  n)2 


(76) 


Subtracting  Equation  64  from  this  last  equation  and  recalling  that 
Ij  =  1  -  0.59n,  the  result  is 


AH 

x 


u  .  h  ,\  H 

H  n  Ap  s  1  \  .  o 

-  +  0  r  — -  I  +  a  — 

2  p  x  /  x 

o  o  /  o 


H  H 

I,  +  A  ~  I,  +  1,  — 
2  o  B  3  5  x 

o 


(77) 


It  is  preferable  to  use  this  equation  to  evaluate  the  average  bottom 
stress  under  the  arrested  saline  wedge. 


PART  X:  EXPERIMENTAL  VALUES  OF  INTERFACIAL  STRESS 


38.  Using  Equations  66  and  77  the  experimental  values  of  average 
interfacial  stress  Tg  and  average  bottom  stress  under  the  wedge  T^ 
are  evaluated  on  the  basis  of  data  shown  in  Tables  1  and  2.  The  evalua¬ 
tion  as  function  of  U  H  /v  are  collected  in  Tables  12  and  13  for  water 

o  o 

depth  Hq  =  23.0  cm  and  45.5  cm,  respectively.  Since  the  values  thus 

arrived  should  be  independent  of  the  ratio  H/B  ,  the  values  of  T  and 

S  2 

Tq  for  the  two  depths  are  shown  in  Figures  11  and  12.  Here,  2Tg/pUo 
is  plotted  logarithmically  against  UqHq/v  and  also  2TQ/pU^  against 
U^H^/v  .  The  scatter  of  the  points  in  both  figures  is  very  large.  This 
is  hardly  surprising  and  was  expected  since  the  values  of  Tg  and  Tq 
are  obtained  in  both  cases  by  subtracting  quantities  of  nearly  equal 
order  of  magnitude,  each  of  these  quantities  separately  being  open  to 
errors  as  large  as  5  percent.  In  particular,  the  fall  of  water  surface 
AH  cannot  be  determined  with  certainty.  In  this  study  AH  is  not 

i 

measured  directly;  it  is  deduced  from  a  curve  of  observed  points  Ah 
versus  x  and  there  is  an  element  of  subjectivity.  Another  cause  for 
the  scatter  is  the  fact  that  the  points  of  the  figures  as  regards  to  n 
are  not  comparable.  By  dimensional  analysis  it  may  be  shown  that 


and 


For  the  data  shown  in  the  figures  n  is  not  constant  but  varies  from 
0.35  to  0.55.  A  good  average  value  of  n  for  all  the  points  is 
n  =  0.45  .  In  Figure  11  a  straight  line  is  drawn  through  the  points, 
with  a  slope  equaling  1.  Half  of  the  observed  points  are  above  the 
line  and  the  other  half  below  the  line;  and  judging  by  the  eye,  the 
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PART  XI:  THEORETICAL  EVALUATION  OF  AVERAGE  STRESSES 


39.  The  evaluation  of  the  stresses  may  be  attempted  on  the  basis 
of  the  measured  velocity  distribution  given  in  Tables  10  and  11. 

It  is  expected  that  this  velocity  pattern  will  apply  almost  to  the  whole 
length  ot  the  arrested  wedge  except  a  stretch  adjacent  to  the  tip. 

40.  Consider  lirst  the  evaluation  of  the  average  interfacial 


stress  T 


The  stress  at  the  interface  is 


du 

T  =  P  , 
s  dz 


z  =  0 


Putting 


and  dividing  by  pU  o 


(J  du/U 

Ts  ^  h  ’  dz/h 
s  s 


du/U 
m  dz/h 


T  /  U  H  Y" 

s  /  o  o  \ 


„  H 
u_  _o 

U  h 
o  s 


Remembering  that  n  =  h  ,/H 

si  o 


x  /U  H  X'1  ,  ,,  h 

s  [  o  o\  1  U  s 

..2  m  l  V  /  n  U  h 
pU  \  /  os 


Introducing  the  average  shear  T 


45 


Placing 


one  now  has 


(86) 


(87) 


Taking  h^/h^  from  Table  4,  Fj  is  computed  by  numerical  integration 
for  a  few  selected  values  of  e  and  n  .  These  values  are  given  in 
Table  14. 

41.  Toward  the  evaluation  of  the  factor  m  one  needs  to  con¬ 
sider  the  velocity  rates  in  the  area  of  the  interface.  The  appropriate 
quantities  are  taken  from  the  entries  of  Tables  10  and  11  and  are  shown 
in  graph  form  in  Figure  13.  The  inclinations  of  the  straight  lines, 
one  for  the  region  of  the  fresh  water  and  the  other  for  the  salt  water, 
are  somewhat  different  in  value;  this  is  in  accordance  with  the  fact  that 
the  viscosity  of  salt  water  is  somewhat  greater  than  that  of  fresh  water. 
On  the  basis  of  these  data  the  factor  m  amounts  to  3.0.  As  the  average 
value'of  n  of  all  the  tests  is  0.45,  the  corresponding  F^  is  read 
from  Table  14  and  for  an  assumed  e  ,  equal  to  0.05,  is  6.53.  Substitu¬ 
tion  of  these  in  Equation  87  yields  for  the  average  interfacial  stress 
the  result 


Figure  13.  Interface  velocity  gradient 
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The  plot  of  this  relation  is  shown  in  Figure  11.  Agreement  between 
the  theory  and  the  observation  is  quite  fair. 

42.  Although  a  slight  mixing  is  present  between  fresh  and  saline 

waters,  this  fact  fails  to  destroy  the  effective  laminar  regime  of  the 

interface.  However,  one  expects  the  initiation  of  a  turbulent  regime 

for  high  Reynolds  number,  but  the  present  tests  fail  to  give  any  indica 

tion  of  the  exact  value  of  U  H  /v  needed.  The  maximum  U  H  /v  noted 

oo  oo 

is  1.2  x  10. 

43.  Flow  in  the  saline  water  close  to  the  bottom  is  turbulent.  < 
the  Reynolds  number  of  flow  is  small,  according  to  Blasius'  hypothesis 
the  velocities  would  vary  as  the  one-seventh  power  of  distance  from 
bottom.  Thus,  one  expects  that 


where  h^  is  the  depth  of  saline  water  and  U  freshwater  current  mean 
velocity  in  a  cross  section.  The  channel  bottom  consists  of  planed  and 
painted  wood.  From  the  experiments  of  Bazin  with  planed  wooden  surfaces 


one  may  deduce  the  velocity  law 


u 

u... 


a  =  5.9 


where  u.,_  =  -^T^/p  •  Dividing  Equation  89  by  Equation  90 


(90) 


U_ 

u_,. 


Solving  for  T  , 
o  ’ 


or 


(91) 


(92) 


Introducing  IM1  =  D(Hq  -  hg)  and  putting 


nH  one  obtains 
o 
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Multiply  the  two  sides  of  this  equation  by  2d £  ,  integrating  between 
£  =  £  and  £  =  1  ,  and  introducing  the  average  bottom  shear  for  the 
entire  length  of  the  wedge  ,  we  now  have 


where 


(93) 


Again  taking  the  values  of  hg/hg^  from  Table  4,  the  definite  integral 
may  be  computed  by  numerical  integration  for  selected  values  of  n  and 
e  .  These  are  given  in  Table  15. 

44.  The  examination  of  lower  velocities  of  the  arrested  saline 
wedge  shown  in  Figure  7  suggests  that  Aq  equals  0.25.  Selecting 
e  =  0.05  and  n  =  0.45  ,  Table  15  gives  =  2.3  .  As  indicated  in 
the  above,  a  =  5.9  .  Inserting  these  in  Equation  73,  the  result  is 


0.018 


(94) 


This  is  plotted  in  Figure  12.  The  estimated  values  of  Tq  are  twice  as 
small  as  the  values  obtained  from  Equation  73.  Probably  the  main  reason 
for  this  large  difference  was  the  obvious  difficulty  of  measuring  small 
velocities  accurately.  In  the  derivation  above  it  was  assumed  tacitly 
that  there  is  similarity  in  velocity  pattern.  In  the  presence  of  mixing 
at  the  interface  this  assumption  would  not  be  valid,  and  the  evaluation 
should  be  carried  out  on  a  different  basis. 
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PART  XII:  DISCUSSIONS 


A  Theory  of  Affine  Shape  of  Arrested  Saline  Wedge 


45.  In  a  study  dealing  with  the  circulation  of  cooling  water 
between  the  intake  and  the  outlet  of  a  thermoelectric  power  plant 
Beta  (1957)  has  shown  that  the  affine  shape  of  an  arrested  wedge  com¬ 
puted  from  the  flow  equations  originally  given  by  Shijf  and  Schoenfeld 
(1953),  is  in  agreement  with  observations.  Comparison  with  theory  is 
made  also  for  observations  from  Keulegan  (1952)  covering  the  range  from 

F  =  0.06  to  F  =  0.40  .  These  fall  on  the  theoretical  curve  corre- 
o  o 

sponding  to  F^  =  0.55  .  Later,  using  a  similar  analysis,  Harleman 
(1961)  has  indicated  that  there  is  agreement  between  theory  and  observa¬ 
tion  as  regards  the  shape  of  arrested  saline  wedges.  These  are  sig¬ 
nificant  findings  and  require  further  consideration. 

46.  The  equation  of  motion  as  relating  to  the  freshwater  layer 

above  the  saline  wedge  may  be  obtained  directly  from  Equation  30. 

Neglecting  the  wall  friction  T  ,  taking  x  equal  to  unity,  writing 

w 

h  for  h  -  h  ,  the  desired  result  is 
w  s 


dh  +  U  dU  _  _  Ts 
dx  g  dx  gph^ 


(95) 


From  Equation  48  neglecting  bottom  friction  Tq  and  placing  P  =  0  (no 
motion  in  the  wedge)  one  has 


A  dh 

dh  A£  _ s 

dx  p  dx 


(96) 


as  relating  to  a  stationary  salt  wedge.  These  agree  with  the  equations 
of  Shijf  and  Schoenfeld  (1953).  Because  of  the  condition  of  continuity 


dh 

U  =  -h  f 
dx  w  dx 


(97) 
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an<1  assuming  that  the  fall  of  surface  water  may  be  neglected 


dh  dh 

w  _  _ s 

dx  dx 


and  then 


dh  ... 

n  s  _  u  dlJ 

U  3 —  =  h  3— 
dx  w  dx 


Introducing  the  latter  in  Equation  95, 


,,  „2  dh 

dh  U  s 

dx  gh  dx 
0  w 


Pgh 


Subtracting  Equation  99  from  Equation  96,  the  result  is 


.  ,,2  \  dh 

Ag  _  U _  j  _ s 

.  P  '  gh  J  dx 


dX  PE  VV  hs/ 


Introducing 


Equation  100  becomes 


A. IT 

1 


Ag  U2  \  dhs  _  Xj  u2  / 1  ,  1  \ 
\P  ’  *Vdx  ’  2  8  Vhw  hs/ 


which  may  be  written  also  as 


P 


U2  \2  H 

o  /  U_\  _o 

"  gH  \  U  /  h 
6  o  \  0/  w 


dh  A.  IT  /,.  \2  / H  H  ' 

_ S  _  _1  O  /  U_  \  j  o  +  o 

dx  "  2  gH  1  U  )  l  h  h 

6  o  \  0/  \  w  s> 


(98) 


(99) 


(100) 


(101  ) 


(102) 
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Plac i ng 
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This  may  be  written  after  placing 


. ,  ,  /.  _2  nL  3  3  4  1  5 

=  V  -  ¥oJ~2  '  \  +  4  nL  "  5  ni 


(107) 


=  r  Fo  r 

o 


(108) 


Here  is  an  expression  which  relates  the  coefficient  of  resistance  of  the 
interface  with  the  length  of  the  arrested  wedge.  It  is  equivalent  to  one 
already  indicated  by  Harleman  (1961).  It  is  valid  only  in  the  case 
where  the  bottom  stress  vanishes  and  the  momentum  of  the  wedge  can  be 
ignored.  Now,  dividing  Equation  105  by  Equation  108  and  writing 
£  =  x/x^  >  as  before,  one  has 


(■  -  3  £  -  ;i  • f  - 


0(nL) 


(109) 


Assuming  that  the  flow  at  the  river  mouth  is  critical 


nr  =  i  - 

L  o 


(110) 


We  note  also  that 


(111) 


Thus,  through  Equations  105,  106,  and  107,  £  may  be  obtained  as  a 


function  of  h  /h  , 
s  si 


for  a  chosen  F 


A  few  determinations  are  shown 


in  Figure  14.  The  circles  represent  experimental  values  for  an  F  of 
0.40  and  are  taken  from  Table  5.  Agreement  between  the  theoretical 
values  and  the  observations  is  good.  In  the  experimental  wedges  friction 
is  present  at  the  flume  walls  and  also  at  the  wedge  bottom.  In  the 
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0  F0  =  0.40 

1  F0  =  0.40 

2  F0  =  0.01 

3  F0  =  0.00 


Figure  14.  Theoretical  affine  shape  of  arrested  wedges 


theoretical  evaluations  these  frictions  are  ignored.  The  observed  close 
agreement  could  be  hardly  expected,  unless  these  two  frictions  are 
of  like  value  and  may  be  imagined  to  be  incorporated  into  the  interfacial 
stress  Tg  .  Another  point  to  remember  is  that  the  critical  flow  con¬ 
dition  at  the  river  mouth,  Equation  110,  is  valid  only  for  values  of  F 
close  to  0.5.  Variation  from  the  expression  becomes  important  as  Fq  ap¬ 
proaches  zero.  In  the  usual  derivation,  it  is  assumed  tacitly  that  at 
the  river  mouth  the  pressure  is  hydrostatic.  When  the  depth  of  salt  wedge 
at  the  river  mouth  is  increased  the  interface  tends  to  be  more  curved, 
affecting  the  pressure  there,  and  this  is  no  longer  hydrostatic.  The 
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Figure  15.  Dependence  of  0  on  F 


critical  conditions  are  not  exactly  known  and  experimental  elucidation 
of  the  matter  is  very  much  desired.  This  was  mentioned  before. 

47.  Equation  108  relates  the  length  of  arrested  saline  wedge  to 
the  friction  coefficient  of  the  interface  A^  .  The  quantity  is 

plotted  against  Fq  in  Figure  15.  For  the  range  between  Fq  =  0.5  and 
F  =0.1  a  good  approximation  for  0(r)r )  is 


(112) 


(|>(Ol)/F“  =  0.015  Fq  ' 

I 

Thus,  a  simple  formula  for  the  resistance  coefficient  of  the  interface 
A^  from  Equation  108  would  be 


10  A.  —  =  0.30  f"5/2 
i  H  o 

o 


(113) 


As  noted  before,  the  length  of  the  arrested  wedges  observed  in  labora¬ 
tory  channels  with  a  depth  width  ratio  H^/B  equal  to  unity  would  be 


or,  since  F  =  V  /V. 

o  r  A 


Introducing  this  in  Equation  113,  practically 

10  A.  =  9.1  F1/2  R'1/2  (114) 

1  o 

Remembering  that  the  mean  of  the  densimetric  Froude  number  of  the  present 
tests  is  Fq  =  0.45  ,  and  inserting  this 

10  A.  =  6.1  (U  H  /v)'1/2  (115) 

i  o  o 

48.  This  would  be  the  value  of  the  coefficient  on  the  basis  of  the 
length  of  wedges  observed  in  the  laboratory  channels  and  the  theoretical 
relation  between  the  coefficient  and  the  length  of  the  wedge,  Equa¬ 
tion  108. 

49.  It  would  be  of  interest  to  compare  the  above  result  with  the 


value  of  A^  to  be  deduced  from  the  observed  in  the  present  study. 

Now,  in  view  of  the  definition  in  Equation  10] 


or 


In  the  tests  of  the  present  study  the  average  value  of  n  is  0.45  and 
the  corresponding  value  of  the  integral  in  the  above  relation  is  1.64. 
Then, 


2T 

— §  =  1.64A 

pU 
^  o 

and  as  the  experiments  of  the  present  study  here  yielded 


the  equation  of  the  upper  line  in  Figure  11,  we  have 


a  result  not  in  agreement  with  Equation  115.  At  the  moment  we  are  unable 
to  discuss  the  reason  for  this  difference. 


Interfacial  Stress  in  Lock  Exchange 

50.  In  a  definitive  study  the  nature  of  interfacial  stress  in  lock 
exchange  was  examined  by  Abraham  and  Eyesink  (1971).  The  coefficient  of 
stress  A^  may  be  defined  as 

A. 

*i  =  " |  PIT  018) 

where  Ur  is  the  relative  velocity  between  the  two  layers  of  the  ex¬ 
change  flow  and  T  the  interface  stress.  Authors  use  the  Weisbach 
coefficient  f  ,  f  =  4A^  .  Resorting  to  an  energy  consideration, 

4A^  is  first  evaluated  utilizing  the  experimental  data  of  Keulegan 
(1957a)  and  Barr  (1963)  on  exchange  flow  involving  laboratory  flumes. 

To  examine  the  effect  of  viscosity  on  the  coefficient  these  authors 
performed  additional  experiments  in  two  channels,  one  small  and  the 
other  3.45  times  as  great  in  the  Delft  Hydraulics  Laboratory.  Evalua¬ 
tion  of  the  coefficient  is  effected  using  the  flow  equations  originally 
formulated  by  Shijf  and  Schoenfeld  (1953).  Initial  conditions  only  are 
considered  coinciding  with  the  densimetric  Froude  number  =  0.9  ,  the 
number  being  defined 


F  = 
r 


(119) 


Hq  being  the  combined  depth  of  the  two  layers.  The  authors  present  the 
evaluations  of  the  coefficient  as  a  function  of  Reynolds  number 


R 

e 


U2R2 


(120) 


where  is  the  current  velocity  in  the  lower  layer  and  R^  is  the 

hydraulic  radius 
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2  2(h  +  B) 

s 


where  hg  is  the  depth  of  the  lower  layer  and  B  the  width  of  channel 


As  the  depth  of  water  Hq  equals  2hg 


■(■A) 


For  the  present  purpose  it  is  desirable  to  express  the  coefficient  as  a 
function  of  the  Reynolds  number  lLhs/v  '  aut^ors'  data  may  be 

readily  expressed  as  a  function  of  U^h  /v  as  well.  This  is  done  in 
Figure  16.  Results  from  the  authors'  experiments  are  in  agreement  with 


O  BARR 
0  KEULEGAN 


Figure  16.  Relation  between  A  and  R  in  exchange  flow 
(after  Abraham  ana  Eyesink  1971) 


the  evaluation  from  the  Barr  and  Keulegan  data.  Three  distinct  regions 
are  discernible.  In  the  laminar  range 


A . =  5 . 5R_ 1  R  <  1000  (123) 

i  e  e 


This  is  also  the  relation  derived  by  Ippen  and  Harleman  (1952)  for  the 
case  of  an  underflow,  that  is,  flow  of  saline  waters  in  an  incline  under 
a  stagnant  pool  of  fresh  water.  In  the  transition  range 


A.  =  0.055R_1/3  1000  <  R  <  40.000  (124) 

l  e  e 

The  authors  point  out  that  the  coefficient  is  independent  of  Reynolds 
number  for  greater  value  of  this  number 

A.  =  0.0017  R  <  40.000  (125) 

l  e 

51.  Agreement  with  our  results  of  the  interfacial  stress  coef¬ 
ficient,  as  expressed  in  Equation  112,  is  lacking.  Now  if  one  would 
accept  for  a  moment  that  the  experimental  procedures  of  our  main  study 
and  the  numerical  analysis  underlying  it  are  admissible,  then  this 
disagreement  could  suggest  that  the  breakdown  of  turbulence  in  the  inter¬ 
faces  of  an  exchange  flow  is  not  similar  to  the  breakdown  in  the  inter¬ 
face  of  an  arrested  saline  wedge. 

Interfacial  Stress  in  Density  Underflow 

52.  A  detailed  study  of  the  mechanism  of  interfacial  flow  between 
fresh  and  saline  waters  was  given  by  .,Loftquist  (  1960).  The  case  con¬ 
sidered  is  the  flow  of  saline  water  under  a  stagnant  body  of  fresh  water 
in  a  horizontal  flume.  The  matters  examined  relate  to  the  interfacial 
gradient  and  stress,  entrainment,  distribution  of  density,  stress  and 
velocity  in  a  cross  section,  and  relation  of  the  length  of  transition 
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layer  to  the  densimetric  Froude  number.  For  the  present,  however,  only 
the  results  of  the  coefficient  of  interfacial  resistance  will  be  con¬ 
sidered.  Loftquist  defines  the  resistance  coefficient  as 


A  = 
m 


T 

m 

2 

pUo 


026) 


where  T  is  in  effect  the  stress  of  interface  and  u  is  the  maximum 
m  o 

velocity  of  saline  waters  in  a  cross  section.  The  densimetric  Froude 
number  is  defined  as 


where  h^  is  hydraulic  radius  of  the  saline  water  layer 


(127) 


bh 

h  -  s 

r  ”  (b  +  h  ) 
s 


(128) 


where  hg  is  the  depth  of  the  saline  water  layer,  2b  the  width,  and 
U  the  mean  velocity  in  the  cross  section  of  the  saline  layer.  In  the 
experiments  conducted,  hg  was  about  18  or  19  cm  and  b  ,  11.5  cm,  so 
that  h^  =  0.37hs  .  The  Reynolds  number  was  defined  as 


R 


L 


Uh 
_ r 

v 

s 


(129) 


where  v  is  the  viscosity  of  saline  waters.  The  stress  I  is  deduced 
s'  m 

from  the  equation  of  motion  of  saline  layer  making  use  of  observed  veloc¬ 
ity  distribution  and  the  fall  of  the  interface.  Evaluations  yield  for 
the  coefficient  of  interfacial  resistance 


A 

m 


4.34R~*P 
L  m 


(130) 
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where  8  is  a  function  of  Rr  and  FT 
m  L  L 

For  the  purpose  of  comparing  with  the  findings  of  the  present  study,  it 

is  desirable  to  express  densimetric  Froude  number  and  Reynolds  number  in 

term  of  h  instead  of  h  ,  thus 
s  r 


U 


k  g  &  h 
P  s 


Uh 

g 

and  R,  =  - 

k  v 


031) 


These  give 


F.  =  0.37F.  ;  R.  =  2.7Rr 
k  L  k  L 


(132) 


Following  the  representation  previously  used,  Equation  97,  the  coef¬ 
ficient  of  interfacial  friction  in  terms  of  mean  velocity  U 


2t 

A.  =  -5 
1  pu2 


Comparing  with  Equation  126 


A,  = 


2  f  \ 


2  2 

During  velocity  traverses  Loftquist  found  uq  =  1.21U  so  that 

A.  =  2.42A 
l  m 


Introducing  in  Equation  130,  A^  from  the  last  relation  and  R^  from 
Equation  132,  the  result  is 


A.  =  28.5RI'18 

i  k  m 


(133) 


The  dependence  of  8  on  Reynolds  number  and  densimetric  Froude  number 

m 

is  given  in  Figure  17.  Construction  of  the  figure  is  original  with  us 
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Figure  17.  Dependence  of  P  on  Reynolds  and  densimetric 
Froude  numbers  (after  Loftquist  1960) 

1/2 

and  is  based  on  data  given  by  Loftquist.  When  F,R,  is  less  than 

K  K 

5,  equals  unity  and  in  this  range  the  interfacial  resistance  coef 

ficient  is 

=  28. 5R'1  (13 

i  k 

a  result  in  close  agreement  with  the  findings  of  the  present  study, 

1/2 

Equation  117.  Further  in  the  range  FR  is  greater  than  4 

P  =  0.25F,R.1/2 
m  k  k 


which  yields 


(135) 


A.  = 
1 


indicating  that  for  larger  values  of  densimetric  Froude  number  the  coef¬ 
ficient  varies  inversely  as  the  one-half  power  of  Reynolds  number.  The 
scatter  of  points  in  Figure  17  renders  those  results  provisional.  In 
fact,  from  Figure  18,  where  is  plotted  against  F^  and  the  scatter 

of  points  is  nearly  of  the  same  extent  as  in  Figure  17,  one  may  also 
conclude  that  the  coefficient  of  interfacial  resistance  varies  directly 
with  F^  and  inversely  with  R^  . 


Figure  18.  Dependence  of  8  on  densimetric  Froude  number 
(after  Loftquist  1960) 
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PART  XIII:  MODIFICATION  OF  STRESS  FORMULA  DUE  TO  ENTRAINMENT 


53.  Equation  66  to  evaluate  the  average  interfacial  stress  T 
was  based  on  the  assumption  that  entrainment  (that  is,  flow  of  saline 
waters  into  the  freshwater  current)  is  absent.  In  the  presence  of 
entrainment  the  coefficients  I  would  undergo  changes.  Possibly  the 
kinetic  reaction  time  will  register  the  greater  change  leading  to  the 
form  I ^  +  Sl^  as  shown  below. 

54.  Denoting  the  entrainments  by  U^  ,  the  vertical  velocity 
components  at  the  free  surface  and  at  the  interface  now  are 


..  dh 

dh  s  ., 

w  =u  —  ,  w .  =  u .  = —  +  U 
s  s  dx  i  i  dx  m 


Then  the  left-hand  term  of  Equation  2 1  changes  to 

b  h  b  h 

.  ,  _  ||  du2  .  /  2  dh  2  dhs 

dz  dy  “  II  d£  dz  dy  »ue  ^  "  11 ; 


J  f(u  IS  +  w  i)  dz  dy  =  JJ 

-b  h  -b  h 


s  dx  i  dx 


U  u.  B 

m  i  j 


and  in  the  place  of  Equation  24  one  now  has 


+b  h 

/  I  (u  ~  +  w  — )  dz  dy  =  j-  (0fU2A)  -  U  u.B 
J  J  \  ax  az )  J  dx  mi 

-b  h 

s 

Integrating  the  right-hand  member  between  x  =  0  and  x  =  x^  and 
dividing  by  B  ,  the  result  is 

x^ 

CfU2h  ,  -  aU2H  -  I  U  u.  dx  (13 

1  wl  oo  I  mi 

o 

where  h  ,  is  the  depth  of  fresh  water  at  river  mouth,  h  ,  =  H  -  h 
wl  r  wl  o  s 

In  view  of  condition  of  continuity 


Hence , 


Assuming  n  =  0.45  ,  adapting  m  =  0.53  ,  assigning  to  U/Uq  the  value 


U 

u 


sufficiently  accurate  for  the  present  purpose  and  carrying  out  the  re¬ 
quired  integration  using  h^/h^  from  Table  4,  we  find 


51,  x 

T7  =  4‘46k  f 

2  o 


(137) 


In  a  previous  investigation  on  mixing  in  arrested  saline  wedges 
(Keulegan  1955a),  it  was  shown  that 

U  =  k(U  -  U  )  k  =  2.12  x  io"4 
m  c 

where  is  the  critical  velocity  for  the  initiation  of  mixing  at  the 

interface  obeying  the  relation 

Uc  =  7.3(vgAp/p)1/3 

In  the  above  analysis  was  neglected  for  simplicity.  Using  k  as 

presently  given 


9.4  x  io'4  f 
o 


(138) 
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PART  XIV:  NEED  FOR  FURTHER  RESEARCH 


56.  Expressions  of  the  affine  shape  of  arrested  saline  wedges 
and  those  of  relative  depth  of  saline  water  at  the  river  mouth  as 
established  by  laboratory  observation  appear  to  be  equally  valid  for 
large  rivers.  As  regards  the  length  of  a  saline  wedge  a  similar  trans¬ 
fer  to  the  prototype  condition  would  not  be  permissible,  especially 
when  the  observations  are  made^ in  narrow  channels.  One  restrictive 
effect  in  such  channels  is  the  friction  of  the  channel  walls.  The  sum¬ 
mary  of  results  from  the  present  and  previous  investigations  was  given 
in  Equation  5  for  H/B  equal  to  1  and  in  Equation  6  for  H/B  equal 
to  2.  To  obtain  the  corresponding  values  for  infinitely  wide  labora¬ 
tory  channels  one  may  resort  to  analysis  using  the  experimentally  de¬ 
termined  values  of  the  interfacial  and  bottom  stresses  (Keulegan  1957b). 
For  a  small  densimetric  Reynolds  of  the  order  of  10  thousand 


(139) 


and  for  large  densimetric  number  of  order  of  10  million 


L 

H 


(HO) 


The  latter  should  also  apply  to  large  rivers  if  the  dependence  of  the  in¬ 
terfacial  and  bottom  stress  on  Reynolds  number  is  the  same  as  in  labora¬ 
tory  small  channels.  This,  however,  is  certainly  open  to  question. 

57.  The  coefficient  of  friction  of  the  interface  obtained  in  this 
study  is  in  agreement  with  Loftquist  data  for  small  velocities  of  saline 
waters  flowing  under  a  pool  of  practically  stagnant  fresh  water  in  a 

horizontal  channel.  When  mixing  is  present,  Loftquist's  data  indicate 

-1/2 

that  the  coefficient  varies  as  R  .  This  effect  is  absent  in  the 
data  of  our  investigation.  The  Delft  Laboratory  tests  on  exchange  flow 
of  lock  operation  reveal  that  during  initial  motion  characterized  by  a 


69 


for 


-1/3 

constant  densimetric  Froude  number  the  coefficient  varies  as  R 
moderately  high  Reynolds  number  and  eventually  at  greater  high  Reynolds 
numbers  above  a  critical  value  it  remains  constant. 

58.  The  direct  determination  of  the  stresses,  interfacial  and 
bottom,  of  an  arrested  wedge  is  a  difficult  matter.  To  reduce  the 
errors  due  to  the  side  frictions  it  is  desired  that  new  investigation 
be  undertaken  with  channels  of  greater  width  in  comparison  with  water 
depths.  Further,  the  channels  should  be  of  such  depth  as  to  allow  flows 
of  larger  Reynolds  number.  In  addition,  closer  attention  should  be  paid 
to  the  matter  of  the  velocities  in  the  area  of  the  wedge  and  in  cross 
sections  more  than  one,  all  uniformly  spaced  across  the  length  of  the 
wedge.  With  the  larger  freshwater  velocities  the  extent  of  mixing  should 
be  determined  and  its  bearing  on  the  wedge  area  velocities  ascertained. 

59.  If  data  are  available  on  the  saline  wedges  in  large  rivers, 
the  question  of  lengths  may  be  readily  examined  by  assuming  that  the 
appropriate  relation  is 


L 

H 
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with  the  constants  A  and  n  to  be  determined.  Taking  the  logarithm 
of  the  two  sides  of  the  equation,  the  resulting  linear  algebraic  equa¬ 
tion  readily  yields  the  values  of  the  unknowns.  The  procedure  certainly 
should  resolve  the  question  if  the  saline  wedge  length  in  rivers  is 
independent  of  Reynolds  number. 
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Table  1 

Arrested  Saline  Wedge  Observed  Data 


H  =  45.5  cm;  B  =  22.9  cm 
o 


No . 

,  cm/sec 

tap  ,  ^ 

cm 

h  ,  ,  cm 
s  1 

0  ,°C 

x  ,  cm 

0 

AH  ,  cm 

IB 

29.3 

0.066 

17.44 

22.8 

2260 

1  .  10 

2B 

25.6 

0.063 

19.36 

22.8 

3230 

1  .  15 

3B 

23.9 

0.068 

21.50 

22.9 

3800 

1.08 

4B 

23.1 

0.066 

20.90 

22.5 

4370 

1.23 

5B 

21.3 

0.058 

22.00 

26.6 

4730 

1  .  15 

6B 

20.5 

0.067 

22.10 

23.0 

5860 

1.22 

8B 

17.9 

0.032 

19.50 

26.7 

2590 

0.53 

9B 

17.0 

0.031 

20.50 

26.0 

2870 

0.54 

10B 

16.3 

0.031 

21.00 

25.3 

3430 

0.56 

1  IB 

15.5 

0.031 

22.50 

24.3 

4170 

0.61 

12B 

14.4 

0.033 

22.50 

25.4 

5040 

0.59 

13B 

14.1 

0.030 

23.0 

25.1 

4970 

0.59 

14B 

13.2 

0.030 

23.75 

25.2 

5440 

0.62 

15B 

14.1 

0.016 

18.00 

26.0 

1680 

0.27 

17B 

12.2 

0.015 

20.5 

26.8 

2768 

0.30 

18B 

11.6 

0.018 

21.0 

26.0 

3290 

0.30 

19B 

11.0 

0.016 

21.75 

26.2 

4680 

0.37 

20B 

10.3 

0.016 

22.7 

25.7 

5260 

0.33 

21B 

9.3 

0.017 

24.0 

26.9 

5590 

0.30 

Table  2 


Arrested  Saline  Wedge  Observed  Data 


H  =  23. 
o 

.  0  cm;  B  = 

22.9  cm 

No. 

U  ,  cm/sec 

a  gm 

cm 

h  ,  ,  cm 
s  1 

6  ,  °C 

x  ,  cm 

0 

AH  ,  cm 

1 

20.  1 

0.0730 

10.0 

23.9 

1510 

0.46 

2 

18.0 

0.0745 

10.0 

24.0 

2776 

0.55 

3 

14.1 

0.0730 

12.0 

24.3 

3809 

0.52 

4 

11.7 

0.0720 

12.7 

24.2 

5632 

0.53 

5 

14.3 

0.0386 

9.9 

24.4 

1130 

0.24 

6 

11.5 

0.0368 

10.9 

24.0 

2172 

0.24 

7 

10.0 

0.0382 

11.1 

22.6 

2624 

0.26 

8 

10.0 

0.0357 

11.3 

23.6 

2577 

0.25 

9 

8.6 

0.0362 

12.6 

23.8 

4117 

0.25 

10 

7.8 

0.0355 

13.0 

23.4 

4957 

0.27 

11 

10.3 

0.0184 

9.7 

24.0 

903 

0.  1 1 

12 

8.2 

0.0181 

10.7 

24.0 

1704 

0.  12 

13 

7.2 

0.0182 

1 1.0 

24.5 

2339 

0.12 

14 

6.2 

0.0178 

11.7 

24.8 

3215 

0.12 

15 

5.1 

0.0172 

12.6 

22.9 

4700 

0.13 

16 

4.8 

0.0174 

13.2 

22.8 

5860 

0.14 

17 

7.2 

0.0092 

9.2 

22.8 

800 

0.055 

18 

5.9 

0.0097 

10.6 

24.5 

1541 

0.050 

19 

4.3 

0.0089 

11.9 

24.8 

2832 

0.058 

20 

3.7 

0.0088 

12.8 

25.9 

3855 

0.057 

21 

3.3 

0.0089 

13.3 

23.7 

5270 

0.065 

I 

f 


Table  5 


1 

1 

Affine 

H  = 
o 

Shape  of  Arrested  Saline  Wedges 

23.0  cm;  B  =  22.9  cm;  F  =  0.40 

0 

Ap/p 

0.0098 

0.0174 

0.0396 

0.0580 

0.0696 

Uq  ,  cm/sec 

5.80 

7.70 

11.9 

14.8 

18.0 

Mean 

t 

h 

s 

/hsl 

0.07 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.05 

0.154 

0.165 

0.137 

0.132 

0.137 

0.145 

0.10 

0.231 

0.231 

0.198 

0.203 

0.187 

0.210 

0.20 

0.330 

0.319 

0.291 

0.286 

0.330 

0.311 

0.30 

0.409 

0.385 

0.368 

0.368 

0.352 

0.376 

0.40 

0.479 

0.434 

0.434 

0.440 

0.418 

0.442 

0.50 

0.547 

0.495 

0.506 

0.506 

0.473 

0.506 

0.60 

0.607 

0.555 

0.577 

0.583 

0.539 

0.572 

0.70 

0.671 

0.621 

0.643 

0.665 

0.605 

0.641 

0.80 

0.742 

0.704 

0.720 

0.748 

0.682 

0.719 

0.90 

0.819 

0.786 

0 . 786 

0.825 

0.875 

0.799 

0.95 

0.886 

0.863 

0.841 

0.885 

0.847 

0.864 

1.00 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

Table  7 

Relative  Surface  Fall  over  Arrested  Salt  Wedge 


H  =  45.5  cm;  B  =  22.9  cm;  F  =0.40 
o  o 


Uq  ,  cm/sec  22.8  14.5  10.6 


Ah/AH 


0, 

.0 

0 

.000 

0. 

.000 

0, 

.000 

0, 

.000 

0. 

.  1 

0 

.051 

0 

.053 

0, 

.  u50 

0, 

.051 

0. 

.2 

0. 

.115 

0, 

.  107 

0. 

,  102 

0. 

.  108 

0. 

.3 

0. 

.  176 

0. 

.173 

0. 

156 

0. 

.  166 

0. 

.4 

0. 

.234 

0. 

,287 

0. 

220 

0. 

.230 

0. 

.5 

0. 

.311 

0. 

,312 

0. 

288 

0. 

,303 

0. 

.6 

0. 

.412 

0. 

.400 

0. 

,368 

0. 

.394 

0. 

.7 

0. 

.508 

0. 

.504 

0. 

.470 

0. 

.494 

0. 

.8 

0 

.621 

0, 

.625 

0. 

.592 

0. 

.612 

0, 

.9 

0 

.770 

0. 

.780 

0. 

.756 

0. 

.766 

1. 

.0 

1. 

.000 

1. 

.000 

1. 

,000 

1. 

,000 

Table  8 

Interface  Velocities 


H  = 
0 

23.0  cm;  B  =  22.9  cm; 

F  =  0.40 

0 

Ap/p 

0.0098 

0.0174 

0.0396 

0.0580 

0.0696 

,  cm/ sec 

5.80 

7.70 

11.9 

14.8 

18.0 

Mean 

C  =  X/X 

U./U 

o 

1 

0.1 

0.425 

0.472 

0.413 

0.468 

0.414 

0.438 

0.2 

0.485 

0.508 

0.471 

0.495 

0.452 

0.482 

0.3 

0.527 

0.520 

0.498 

0.509 

0.484 

0.507 

0.4 

0.555 

0.534 

0.519 

0.515 

0.504 

0.525 

0.5 

0.574 

0.541 

0.523 

0.520 

0.516 

0.535 

0.6 

0.594 

0.546 

0.524 

0.523 

0.522 

0.542 

0.7 

0.594 

0.551 

0.522 

0.510 

0.518 

0.539 

0.8 

0.602 

0.550 

0.515 

0.499 

0.511 

0.534 

0.9 

0.609 

0.552 

0.500 

0.485 

0.500 

0.529 

1.0 

0.590 

0.550 

0.475 

0.470 

0.435 

0.502 

Table  9 

Surface  Velocities 


H  =  23.0  cm;  B  =  22.9  cm;  F  =  0.40 
o  o 


Ap/p 

0.0098 

0.0174 

0.0396 

0.0580 

0.0696 

Mean 

,  cm/sec 

5.80 

7.70 

22.9 

14.8 

18.0 

C  =  x/xo 

u  /u 

s 

0.0 

1.032 

1.043 

1.024 

1.055 

1 .048 

1.040 

0.1 

1.042 

1.051 

1.023 

1.060 

1.055 

1.046 

0.2 

1.051 

1.057 

1.030 

1.070 

1.062 

1.054 

0.3 

1.063 

1.065 

1.034 

1.080 

1.060 

1.060 

0.4 

1.072 

1.069 

1.042 

1.080 

1.060 

1.065 

0.5 

1.077 

1.076 

1.050 

1.075 

1.060 

1.066 

0.6 

1.081 

1.078 

1.058 

1.075 

1.058 

1.070 

0.7 

1.084 

1.083 

1.066 

1.070 

1.048 

1.070 

0.8 

1.086 

1.086 

1.073 

1.060 

1.035 

1.068 

0.9 

1.089 

1.090 

1.071 

1.045 

1.020 

1.063 

1.0 

1.091 

1.095 

1.070 

1.035 

1.000 

1.058 

Table  10 


Velocity  Distribution  in  Salt  Wedge 

H  =  23.0  cm;  B  =  22.9  cm;  F  =0.40 
o  o 


Series 

1 

2 

3 

4 

5 

Mean 

Ap/p 

0.0098 

0.0174 

0.0396 

0.0580 

0.0696 

h  ,  cm 
w 

15.1 

14.9 

15.6 

16.4 

17.4 

15.8 

h 

s 

7.4 

7.7 

7.1 

6.7 

6.2 

7.0 

z  ’  /  h  * 
s 

u/U 

0.00 

0.540 

0.536 

0.567 

0.480 

0.527 

0.530 

-0.04 

0.429 

0.406 

0.438 

0.383 

0.406 

0.412 

-0.08 

0.330 

0.358 

0.351 

0.312 

0.328 

0.336 

-0.12 

0.257 

0.216 

0.287 

0.258 

0.245 

0.253 

-0.16 

0.202 

0.163 

0.235 

0.212 

0.184 

0.199 

-0.20 

0.152 

0.129 

0.  195 

0.171 

0.140 

0.157 

-0.28 

0.060 

0.066 

0.135 

0.109 

0.074 

0.089 

-0.36 

-0.040 

-0.041 

0.072 

0.070 

0.027 

0.018 

-0.44 

-0.108 

-0.086 

-0.015 

-0.007 

-0.042 

-0.052 

-0.52 

-0.148 

-0.129 

-0.079 

-0.081 

-0.117 

-0.111 

-0 . 64 

-0.182 

-0.157 

-0.135 

-0.141 

-0.155 

-0. 154 

-0.68 

-0. 189 

-0.161 

-0.147 

-0.152 

-0.160 

-0.162 

-0.72 

-0.192 

-0.162 

-0.153 

-0.156 

-0.162 

-0.165 

o 

1 

-0.194 

-0. 161 

-0.154 

-0. 156 

-0. 160 

-0.165 

-0.80 

-0.194 

-0.159 

-0.152 

-0.152 

-0. 155 

-0.162 

-0.84 

-0.191 

-0.155 

-0. 148 

-0.143 

-0.145 

-0.156 

-0.90 

-0. 177 

-0. 140 

-0. 132 

-0.117 

-0. 1 12 

-0.136 

-0.96 

-0.138 

-0.103 

-0.095 

-0.071 

-0.062 

-0.094 

-1.00 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

z 


is  measured  from  interface  upward. 


Table  11 

Velocity  Distribution 

in  Fresh  Water  over  the  Salt  Wedge 

H  =  23.0  cm;  B  =  22.9  cm;  F  =  0.40 
o  ’  o 


Series 

1 

2 

3 

4 

5 

Mean 

Ap/p 

0.0098 

0.0174 

0.0396 

0.0580 

0.0696 

h  ,  cm 
w 

15.1 

14.9 

15.6 

16.4 

17.4 

15.8 

h 

s 

7.4 

7.7 

7.1 

6.7 

6.2 

7.0 

z’/h* 

u/U 

0.00 

0.540 

0.536 

0.567 

0.480 

0.527 

0.530 

0.04 

0.648 

0.666 

0.665 

0.580 

0 . 646 

0.641 

0.08 

0.736 

0.746 

0.735 

0.651 

0.708 

0.715 

0.12 

0.801 

0.792 

0.785 

0.706 

0.755 

0.768 

0.16 

0.846 

0.827 

0.822 

0.746 

0.789 

0.806 

0.20 

0.880 

0.857 

0.852 

0.780 

0.830 

0.840 

0.24 

0.908 

0.883 

0.876 

0.810 

0.852 

0.866 

0.32 

0.950 

0.921 

0.915 

0.855 

0.886 

0.905 

0.40 

0.982 

0.952 

0.944 

0.889 

0.922 

0.938 

0.56 

1.023 

0.997 

0.987 

0.956 

0.966 

0.986 

0.72 

1.045 

1.026 

1.015 

1.018 

0.994 

1.020 

1.04 

1.062 

1.064 

1.052 

1.056 

1.036 

1.054 

1.20 

1.060 

1.070 

1.064 

1.064 

1.050 

1.062 

1.36 

1.052 

1.072 

1.071 

1.074 

1.062 

1.066 

1.52 

1.040 

1.077 

1.070 

1.080 

1.074 

1.068 

1.68 

1.038 

1.058 

1.064 

1.080 

1.078 

1.060 

1.84 

1 . 026 

1.046 

1.053 

1.078 

1.078 

1.056 

2.00 

1.011 

1.038 

1.038 

1.072 

1.070 

1.047 

2.32 

0.991 

1.038 

1.021 

1.072 

1.063 

1.037 

* 


z'  is  measured  from  interface  upward. 


Table  12 

Evaluation  of  Total  Interfacial  and  Bottom 

Frictional  Stress  Forces 

H  =  23.0  cm:  B  =  22.9  cm;  v  =  0.00919 
o 


No. 

Ap/p 

2V  /V. 
r  A 

h  /H 
s  0 

U  H  /v 

0  o 

2T  /pU2 
s  r  0 

2T  /pU2 
0  0 

1 

0.0730 

1.001 

0.370 

49.2  x  103 

10  x 

io"4 

33  x  io 

2 

0.0745 

0.891 

0.405 

45.2 

7 

31 

3 

0.0730 

0.698 

0.480 

35.4 

16 

22 

4 

0.0720 

0.588 

0.523 

16.8 

20 

11 

5 

0.0386 

0.972 

0.380 

35.9  x  io3 

19  x 

10~4 

30  x  io 

6 

0.0368 

0.806 

0.425 

28.8 

43 

32 

7 

0.0382 

0.682 

0.485 

25.1 

56 

13 

8 

0.0357 

0.710 

0.475 

25.1 

44 

11 

9 

0.0362 

0.601 

0.530 

21.6 

27 

27 

10 

0.0355 

0.542 

0.540 

18.5 

48 

15 

11 

0.0184 

1.015 

0.372 

25.8  x  103 

20  x 

io"4 

43  x  io 

12 

0.0181 

0.773 

0.450 

21.8 

12 

68 

13 

0.0182 

0.756 

0.455 

18.0 

48 

13 

14 

0.0172 

0.616 

0.510 

15.4 

27 

52 

15 

0.0172 

0.523 

0.545 

12.9 

130 

26 

16 

0.0174 

0.479 

0.570 

11.9 

38 

22 

17 

0.0092 

1.001 

0.370 

18.1  x  io3 

56  x 

io'4 

50  x  15' 

18 

0.0097 

0.803 

0.435 

14.7 

56 

41 

19 

0.0089 

0.616 

0.510 

10.7 

56 

49 

20 

0.0088 

0.532 

0.510 

9.3 

80 

23 

21 

0.0089 

0.472 

0.570 

8.3 

32 

24 
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Table  13 

Evaluation  of  Total  Interfacial  and  Bottom 

Frictional  Stress  Forces 

H  =  45.5  cm;  B  =  22.9  cm;  v  =  0.00950 
o 


No . 

Ap/p 

2V  /V. 
r  A 

h  ,/H 
si  o 

U  H  /v 

O  0 

2Veuo 

2T  /pin 

0  *0 

4 

-4 

-4 

1 

0.066 

1.08 

0.382 

14.3  x  10 

4  x  10 

20  x  10 

2 

0.063 

0.89 

0.425 

13.7 

6 

40 

3 

0.068 

0.86 

0.472 

11.4 

6 

52 

4 

0.066 

0.85 

0.459 

11.0 

4 

20 

5 

0.058 

0.84 

0.483 

10.2 

2 

18 

6 

0.067 

0.76 

0.485 

9.8 

6 

24 

,  4 

-4 

-4 

8 

0.032 

0.95 

0.427 

8.6  x  io 

4  x  io 

42  x  10 

9 

0.031 

0.91 

0.450 

8.2 

4 

36 

10 

0.031 

0.89 

0.461 

7.8 

2 

28 

11 

0.031 

0.87 

0.493 

7.4 

8 

32 

12 

0.033 

0.80 

0.493 

6.9 

10 

32 

13 

0.030 

0.75 

0.502 

6.7 

16 

18 

14 

0.030 

0.76 

0.521 

6.3 

16 

16 

4 

-4 

-4 

15 

0.016 

1.03 

0.395 

6.7  x  10 

4  x  io 

o 

X 

o 

17 

0.015 

0.92 

0.450 

5.8 

16 

16 

18 

0.018 

0.80 

0.461 

5.5 

20 

28 

19 

0.016 

0.80 

0.477 

5.3 

8 

16 

20 

0.016 

0.75 

0.498 

4.9 

12 

28 

21 

0.017 

0.65 

0.526 

4.4 

12 

60 
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Table  14 

Values  of  the  Definite  Integral 


n 

£ 

0.8 

0.7 

0.6 

F1 

0.5 

0.4 

0.3 

0.05 

4.74 

4.94 

5.33 

5.99 

7.07 

8.88 

0.10 

4.26 

4.19 

4.70 

5.25 

6.15 

7.71 

0.15 

3.91 

3.80 

4.25 

4.73 

5.51 

6.88 

0.20 

3.60 

3.46 

3.87 

4.28 

4.97 

6.19 

APPENDIX  A:  NOTATION 


Half  width  of  channel 
Width  of  channel 

Densimetric  Froude  number  F,  =  gh  for  saline 

underflow 

Densimetric  Froude  number  F,.  =  gh  for  saline  under- 

T  *  0  ^  IT 

flow  below  stagnant  fresh  water  H 

Densimetric  Froude  number  F  =  gH  arrested  saline 

2  -  o  o  p  o 

wedge.  F  =  U  /  — gH  for  lock  exchange  flow 
o  r  p  °  o  ® 

Total  depth  of  water  in  the  area  of  the  arrested  wedge 

Hydraulic  radius  of  saline  waters  of  a  saline  underflow, 

h  =  bh  / (b  +  h  ) 
r  s  s 

Depth  of  saline  water  in  an  arrested  saline  wedge;  depth  of 
saline  water  of  an  underflow  or  the  depth  of  saline  water  in  a 
lock  exchange  flow 

Depth  of  saline  water  of  an  arrested  wedge  at  the  river  mouth 
Depth  of  fresh  water  in  the  layer  over  the  arrested  wedge 
Depth  of  fresh  water  at  the  river  mouth 


Depth  of  fresh  water  at  the  tip  of  an  arrested  saline  wedge; 
total  depth  of  the  fresh  and  the  saline  water  layers  of  a  lock 
exchange  flow 

Total  length  of  arrested  saline  wedge.  Same  as  x 


M  Momentum  of  body  liquid  of  arrested  saline  wedge,  and  of  body 
of  fresh  water  resting  on  the  entire  length  of  arrested  saline 
wedge 

+  M  Value  of  M  at  time  t 

o 

+  M  Value  of  M  at  time  t  +  At 

o 

n  Ratio  h  /h  , 

s  si 

P.  Pressure  total  force  acting  on  the  entire  interface 

Pj.P^  Pressure  forces  on  the  upstream  and  the  downstream  faces 
in  freshwater  layer 

Q  Discharge  of  river  waters 

R  Reynolds  number  U  H  /v  for  arrested  wedge 


A1 


R^  Reynolds  number  ^2^2^V  ^or  exc^an8e  flow 

R^  Reynolds  number  Uh^/v  for  saline  underflow 

R^  Reynolds  number  Uh^/v  for  saline  underflow 

R^  Hydraulic  radius  of  lock  exchange  flow,  hsB/2(hg  +  B) 

T  Mean  bottom  stress  averaged  along  the  entire  length  of  the 

arrested  saline  wedge 

T  Mean  interfacial  stress  averaged  along  the  entire  length  of 
the  arrested  saline  wedge 

T  Retarding  force  from  the  sidewalls  in  the  entire  freshwater 

w  layer 

T  Retarding  force  from  the  sidewalls  in  the  entire  saltwater 

w2  , 

layer 

u  Particle  velocity;  root  mean  square  of  velocities  in  strip 
plane 

u  Maximum  velocity  in  the  saline  layer  of  an  underflow 

o  * 

U  Mean  velocity  of  fresh  water  in  the  layer  over  an  arrested 
saline  wedge;  mean  velocity  of  saline  water  in  the  layer 
under  the  stagnant  water  of  an  underflow 

IK  Velocity  of  points  at  the  interface 

Um  Velocity  of  transport  across  the  interface  due  to  mixing 

U  Mean  velocity  of  fresh  water  at  the  tip  of  an  arrested  saline 
wedge.  Same  as  V^ 

Relative  velocity  of  layers  in  a  lock  exchange  flow 
Velocity  of  surface  waters  in  an  arrested  saline  wedge 
Uj  Mean  velocity  of  fresh  water  at  the  river  mouth 

Mean  velocity  of  saline  layer  in  a  lock  exchange  flow 

V  River  velocity  (V  =  U  ) 

r  J  v  o 

V^  Densimetric  velocity,  V^  -  gHc 

x  Coordinate  in  the  direction  of  freshwater  flow,  measured  from 
the  tip  of  an  arrested  saline  wedge 

xq  Length  of  an  arrested  saline  wedge.  Same  as  L 

y  Coordinate  in  the  channel  transverse  direction,  measured  from 
channel  midvertical  plane 

z,z'  Coordinates  in  the  vertical;  z  measured  from  the  bottom  of 
horizontal  channel;  z'  from  the  interface 

a  Boussinesq  coefficient  of  velocity  distribution  in  relation  to 
freshwater  flow  over  arrested  saline  wedge 

A2 


Pm 

Ah 

AH 


A. 

1 


o 

M 


T  . 
1 


Boussinesq  coefficient  in  relation  to  saline  waters  of  arrested 
wedge 

Numerical  factor  in  Loftquist  formula  for  A 

II: 

Fall  of  free  surface  in  the  area  of  an  arrested  saline  wedge 

Fall  of  free  surface  at  the  river  mouth,  H  =  AH  +  h  ,  +  h  .  ) 

o  wl  si 

Relative  depth  of  salt  wedge,  h^/HQ 

Relative  depth  of  salt  wedge  at  river  mouth,  same  as 

0  =  h  ,/H 
s  1  o 

Temperature  of  water  expressed  in  centigrade 
A  function  of  r|  defined  by  Equation  107 
Coefficient  of  resistance  of  river  channel 

2 

Interfacial  coefficient  of  resistance.  A.  =  2t  /pU  for 

1  1  2 

arrested  saline  wedges  and  for  saline  underflow.  A.  =  2l./pU 
for  lock  exchange  flow 

2 

Interfacial  coefficient  of  resistance  A  =  T  /pU  for  saline 
.  m  m  r  o 

underflow 

Coefficient  of  resistance  defined  by  Equation  65 
Viscosity  of  water 
Kinematic  viscosity  of  water 
Kinematic  viscosity  of  saline  water 

Interfacial  stress  for  arrested  saline  wedge  or  for  lock  ex¬ 
change  flow  or  for  underflow 

Maximum  value  of  stress  in  underflow,  practically  same  as 
Bottom  stress  of  an  arrested  saline  wedge 
Interfacial  stress 


In  accordance  with  letter  from  DAEN-RDC,  DAEN-ASI  dated 
22  July  1977,  Subject:  Facsimile  Catalog  Cards  for 
Laboratory  Technical  Publications,  a  facsimile  catalog 
card  in  Library  of  Congress  MARC  format  is  reproduced 
below. 
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